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ABSTRACT 

The hypothesis that ultra-high energy (> 10 19 eV) cosmic rays (UHECRs) are accel- 
erated by gamma-ray burst (GRB) blast waves is assumed to be correct. Implications 
of this assumption are then derived for the external shock model of gamma-ray bursts. 
The evolving synchrotron radiation spectrum in GRB blast waves provides target pho- 
tons for the photomeson production of neutrinos and neutrons. Decay characteristics 
and radiative efficiencies of the neutral particles that escape from the blast wave are 
calculated. The diffuse high-energy GRB neutrino background and the distribution 
of high-energy GRB neutrino events are calculated for specific parameter sets, and a 
scaling relation for the photomeson production efficiency in surroundings with different 
densities is derived. 

GRBs provide an intense flux of high-energy neutrons, with neutron-production ef- 
ficiencies exceeding ~ 1% of the total energy release. The radiative characteristics of 
the neutron /?-decay electrons from the GRB "neutron bomb" are solved in a special 
case. Galaxies with GRB activity should be surrounded by radiation halos of ~ 100 
kpc extent from the outflowing neutrons, consisting of a nonthermal optical/X-ray syn- 
chrotron component and a high-energy gamma-ray component from Compton-scattered 
microwave background radiation. The peak luminosity emitted by the diffuse /3-electron 
halo from a single GRB with > 2 x 10 53 ergs isotropic energy release is ~ 10 35 ergs 
s j with a potentially much brighter signal from the neutron-decay protons. The de- 
cay halo from a single GRB can persist for > 0.1-1 Myr. Stronger neutrino fluxes and 
neutron decay halos can be produced by external shocks in clumpy external media and 
in scenarios involving internal shock scenarios, so detection of neutrinos associated with 
smooth-profile GRBs could rule out an implusive GRB central engine and an external 
shock model for the prompt phase. 

The luminosity of sources of GRBs and relativistic outflows in L* galaxies such as 
the Milky Way is at the level of ~ 10 40±1 ergs s _1 . This is sufficient to account for 
UHECR generation by GRBs. We briefly speculate on the possibility that hadronic 
cosmic rays originate from the subset of supernovae that collapse to form relativistic 
outflows and GRBs. 
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Subject headings: cosmic rays — galaxies: halos — gamma rays: theory — gamma rays: 
burst 

1. Introduction 

The distance scale to the sources of gamma-ray bursts (GRBs) with durations > 1 s has been 
established as a consequence of observations made with the Beppo-SAX satellite (Costa et al. 
1997; van Paradijs et al. 1997). The Beppo-SAX discovery of decaying X-ray afterglows permits 
follow-on optical observations that give redshift determinations from absorption and emission lines 
in optical transient counterparts or from directionally coincident host galaxies. Nearly 20 GRB 
sources have measured redshifts (for a recent review, see van Paradijs et al. (2000)), with a mean 
redshift z ~ 1 for the sample. The distribution of redshifts is as yet poorly established, but ranges 
from z = 0.0085 for GRB 980425 to z = 4.50 for GRB 000131. The redshift of GRB 980425 is based 
upon its temporal and spatial coincidence with SN 1998bw (Galama et al. 1998; Kulkarni et al. 1998; 
Pian 2000), and points to a relationship between GRB sources and supernovae (SNe). A GRB/SN 
relationship is strengthened by the detection of highly reddened excesses in the optical afterglows 
of several GRBs, which would arise if supernova ejecta, powered by the decay of radioactive 56 Ni 
(Bloom et al. 1999; Reichart 1999; Galama et al. 2000), are formed in GRB explosions. Measured 
apparent isotropic 7-ray energy releases range from E 1 ~ 10 48 ergs for GRB 980425 to ~ 2.4 x 10 54 
ergs for GRB 990123 at z = 1.60, with £ 7 >3x 10 51 ergs in all cases except GRB 980425 (Frail 
et al. 2001). Achromatic temporal breaks in the optical light curves of GRB 990123 (Kulkarni et 
al. 1999) and GRB 990510 (Harrison et al. 1999) suggest, however, that the most luminous GRBs 
might be beamed, so that only directional energy releases are actually measured. In the case of 
GRB 990123 (Briggs et al. 1999), the directional 7-ray power and 7-ray energy release reach peak 
values dLj/dQ, ~ 3 x 10 51 ergs s _1 sr _1 and dE/dVL ~ 2 x 10 53 ergs sr _1 , respectively. 

Considerable evidence linking the sources of GRBs with star-forming regions in galaxies has 
recently been obtained (e.g., Lamb 1999; Djorgovski et al. 2001). Optical transients associated 
with GRBs are superposed on the stellar fields of associated host galaxies in essentially all 14 
cases of GRBs with deep follow-up optical observations (van Paradijs et al. 2000; Fruchter et al. 
1999; Bloom et al. 1999a; Odewahn et al. 1998), rather than far outside the galaxies' disks, as 
might be expected in a scenario of merging neutron stars and black holes (Narayan et al. 1992). 
Host galaxies that are directionally coincident with optical transients discovered within the field 
of GRB X-ray afterglows have blue colors, consistent with galaxy types that are undergoing active 
star formation (Fruchter et al. 1999; Castander and Lamb 1999a,b). The host galaxy luminosities 
are consistent with a Schechter luminosity function (Schaefer 2000), and span a wide range of 
extinction-corrected R magnitudes from R ~ 13 for the host galaxy of GRB 980425 associated with 
SN1998bw to R > 27.1 for GRB 980326 (Schaefer 2000; Hogg and Fruchter 1999). Lack of optical 
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counterparts in some GRBs such as GRB 970828 and GRB 991226, which have associated radio 
counterparts (Frail et al. 1999), could be due to extreme reddening from large quantities of gas and 
dust in the host galaxy (e.g., Owens et al. 1998). X-ray evidence (Piro et al. 2000) for Fe K a -hne 
signatures in GRB 991216, requiring large masses and column densities of nearby gas (Bottcher 
2000), also indicates that GRBs originate in regions with active star formation. 

Knowledge of the distance scale to GRBs makes it possible to determine their effects on the 
surrounding environment. Some of the claimed effects of GRB explosions are the formation of HI 
shells and stellar arcs (Efremov et al. 1998; Loeb and Perna 1998), the melting of dust grains by 
GRB UV radiation to produce flash-heated chondrules in the early Solar system (McBreen and 
Hanlon 1999), and the formation of sites of enhanced annihilation radiation in the interstellar 
medium (ISM) originating from large numbers of mildly relativistic positrons produced by a GRB 
(Dermer and Bottcher 2000; Furlanetto and Loeb 2002). UV and X-rays from nearby GRBs could 
also have produced biologically significant dosages on Earth in the past (Scalo and Wheeler 2002). 

Another effect of GRBs, proposed prior to the Beppo-SAX discovery, is that GRB sources 
accelerate the highest energy cosmic rays. Milgrom and Usov (1995) argued for this connection on 
the basis of a directional association of two > 10 20 eV air shower events with earlier BATSE GRBs. 
Waxman and Coppi (1996) pointed out, however, that the intergalactic field must disperse the 
arrival time of the cosmic rays by > 50 yrs to be consistent with the detection rate of GRBs. Vietri 
(1995) noted that the isotropy of the UHECR arrival direction was consistent with the isotropic 
distribution of GRB sources, and that the extreme energies of UHECRs could be explained through 
first-order Fermi acceleration by a relativistic blast wave with Lorentz factor T. At each shock 
crossing, a particle would increase its energy by a factor ~ 4T 2 ~ 4 x 10 5 (r/300) 2 , so that only a 
few such cycles would suffice to produce UHECRs starting from low-energy particles. The efficiency 
to accelerate low-energy particles to ultra-high energies through relativistic shock acceleration has 
since been shown to be infeasible (Gallant and Achterberg 1999; Gallant et al. 1999). Following the 
first shock crossing, the blast wave intercepts the particle before its angular deflection from the shock 
normal is much larger than 1/T; thus subsequent cycles lead to energy increases by only factors 
of ~ 2. Second-order Fermi acceleration, for example, due to magnetohydrodynamic turbulence 
generated by charged dust or irregularities in the external medium (Waxman 1995; Schlickeiser 
and Dermer 2000; Dermer and Humi 2001), or by first-order Fermi acceleration involving putative 
shocks in a relativistic wind (Waxman 1995) could, however, accelerate UHECRs in GRB blast 
waves. 

Both Vietri (1995) and Waxman (1995) pointed out a remarkable coincidence between the 
energy density of the highest energy cosmic rays and the power of GRB sources within the Greisen- 
Zatsepin-Kuzmin (GZK) radius, outside of which UHECRs are degraded by photomeson production 
on the cosmic microwave background. If GRB sources convert a comparable amount of energy into 
UHECRs as is detected in the form of 7 rays, then these sources can account for the observed 
intensity of UHECRs. The comparisons of Vietri (1995) and Waxman (1995) made use of statistical 
studies where the most distant GRBs detected with BATSE were assumed to be at z ~ 1. Redshift 
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measurements of GRB sources now permit more refined studies of GRB statistics, yielding the 
comoving space density of GRB sources and the volume-averaged energy injection rate of GRB 
sources into the ISM. This coincidence can therefore be more carefully tested. 

In this paper, it is assumed that the sources of UHECRs are GRBs. We then examine the 
implications that follow from this assumption. Theoretical problems with accelerating particles to 
ultra-high energies are not dealt with here (see, e.g., Rachen and Mesazaros (1998); Vietri (1998b); 
Dermer (2001); Dermer and Humi (2001)). In Section 2, we summarize a recent statistical study 
employing the external shock model for GRBs (Bottcher and Dermer 2000) and compare it with 
other statistical studies of GRBs and the constant-energy-reservoir result of Frail et al. (2001). 
An external shock model is more energetically efficient than internal shocks to generate 7 rays 
in the prompt phase of a GRB, so this study yields a lower limit to the energy production rate 
of GRB sources per comoving volume. We show that an external shock model is consistent with 
the UHECR/GRB hypothesis, so that the coincidence originally identified by Vietri (1995) and 
Waxman (1995) holds. 

In Section 3, the evolving temporal and spectral behavior of synchrotron radiation in GRB 
blast waves is characterized. This radiation provides a target photon source for high-energy pro- 
tons, and we calculate neutron and neutrino production from photopion processes in GRB blast 
waves. Neutral particle production spectra, integrated over the prompt and afterglow phases of a 
GRB, are calculated. The diffuse high-energy neutrino background and the distribution of neu- 
trino event rates are calculated in Section 4. The outflowing neutrons decay to form high-energy 
protons and electrons. In Section 5, the radiation halos formed through synchrotron and Thomson 
processes of neutron /3-decay electrons are derived in the special case of a power-law distribution of 
neutrons that are impulsively released from a GRB source. The hypothesis that hadronic cosmic 
rays originate from the subset of supernovae that collapse to form relativistic outflows and GRBs is 
briefly considered in Section 6. Fuller discussions of this hypothesis can be found elsewhere (Der- 
mer 2000a, b). Summary and conclusions are given in Section 7. Appendix A gives the synchrotron 
radiation limit used to determine the maximum proton energies, and Appendix B gives a scaling 
for the photomeson production efficiency in surrounding circumburster media (CBM) with different 
densities. 

2. Statistics and Energetics of UHECRs and GRBs 

2.1. GRB Statistics 

The BATSE instrument on the Compton Gamma Ray Observatory provides a data base of 
peak count rates and peak fluxes for several thousand GRBs with unknown redshifts (Paciesas et 
al. 1999). Many attempts have been made to derive the GRB rate density and mean luminosities by 
modeling this size distribution. Even constraining the implied redshift distribution to be consistent 
with the z-distribution for the GRBs with measured redshifts, it has not been possible to derive 
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these quantities unambiguously from the size distribution alone. Uncertainties in determining the 
rate density of GRBs arise from lack of knowledge of the redshift distribution (Totani 1997), the 
luminosity function (Mao and Mo 1998; Krumholz et al. 1998; Hogg and Fruchter 1999; Schmidt 
1999), and the spectral shape (Mallozzi et al. 1996) of GRBs. A useful simplification (Totani 1997; 
Wijers et al. 1998; Totani 1999) is to assume that the GRB rate density is proportional to the 
star formation rate (SFR) history of the universe as traced, for example, by faint galaxy data in 
the Hubble Deep Field (Madau et al. 1998) (which may, however, seriously underestimate the true 
star formation rate at z > 1 (Blain et al. 1999)). An important result is that GRBs are unlikely 
to be standard candles, whether or not their birth rate follows the SFR or a range of reasonable 
evolutionary models (Schmidt 1999; Hogg and Fruchter 1999; Krumholz et al. 1998). 

To constrain the models further, Bottcher and Dermer (2000) jointly modeled the distributions 
of peak flux, duration, and peak photon energies of the vF v spectra of GRBs using an analytic 
representation (Dermer et al. 1999) of temporally evolving GRB spectra in the external shock 
model of GRBs. The assumption that the GRB source density followed the star formation history 
of the universe was maintained, and a flat ACDM cosmology with (f^o, ^a) = (0.3, 0.7) and Hubble 
constant Hq = lOO/i km s _1 Mpc -1 , with h = 0.65, was used. The model flux was folded through 
the simulated triggering response of a BATSE detector to determine detect ability. This approach 
requires that the total energy Eq and the initial blast wave Lorentz factor Tq of a GRB source be 
specified. The burst luminosity is then calculated through the standard blast-wave physics that 
yielded the analytic representation of the GRB spectrum. The analytic model is degenerate in the 
quantity uqTq, where no is the density of the surrounding medium, which is assumed to be uniform. 
The photomeson production efficiency can be scaled from no, as shown in Appendix B. 

Bottcher and Dermer (2000) showed that fixed values of Eq and Tq could not explain the 
observed distributions, and that broad ranges of values are required. The comoving differential 
density distribution of GRB sources was obtained by assuming that the Eq and Tq distributions 
are separable from the redshift distribution and are adequately described by single power-law 
distributions. The rate-density distribution nGRB^Cb To! z) of GRB sources that gives a reasonable 
fit to the size, duration, and peak photon energy distributions, in units of Gpc~ 3 yr _1 E^ Tq 1 , is 

n GRB (E 52 , Tq; z) = 0.022 E(z) ^^Fo - 25 H[E 52 ; HT 4 , 10 2 ] H[T ; 1, 260] . (1) 

In equation (1), Eq = W 52 E^2 ergs, and the Heaviside function is defined such that H[x;a,b] = 1 
for a < x < b, and H[x;a, b] = otherwise. The range of To given here corresponds to a density 
no = 10 2 cm -3 although, again, the model is degenerate in the quantity noFg. The analytic 
representation of the SFR function, normalized to unity at z = 0, is 

' 1 , for z < 0.3 

, 5 • 10 2 - 1 , for 0.3 < z < 1.1 

E < 2 >= 6.3, for 1.1 < 3 i 2.8 (2) 

, 210 • l0-°- 4 ( 2+1 ) , for 2.8 < z < 10. 
(Note that the z < 0.3 branch of this function was omitted in equation (12) of Bottcher and Dermer 
(2000).) 
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The burst rate and energy release rate per unit comoving volume by GRB progenitors can be 
easily obtained from equation (1). In the local universe, we find that 



n G RB(z = 0) = / dE 52 / dT h GRB (E 52 ,T ; z = 0) * 3.6(£ 5 fT°- 52 Gpc~ 3 yr" 1 (3) 



for the burst rate density, where Egg 71 is the minimum energy of GRB sources in units of 10 52 
ergs, and the expression on the right-hand-side of this equation is valid when E^ n <C 100. The fit 
to the BATSE statistics is not sensitive to the value of E™ in when E™ in <C 1. When E™ in = 1, 
h GRB (z = 0) ^ 3.6 Gpc~ 3 yr -1 . If GRB 980425 is assumed to be associated with SN 1998bw, then 
the statistical model requires that E^ n = 10 -4 . In this case, tigrb{z = 0) = 430 Gpc~ 3 yr -1 , and 
most GRBs have low energy and luminosity and are consequently not observed. The event rate is 
therefore very sensitive to the number of faint bursts which is not well-constrained by present data 
(Mao and Mo 1998). Beaming will increase the rate density of sources by the inverse of the mean 
beaming fraction compared to the isotropic value given here. 

The local energy emissivity of the sources of GRBs, from equation (1), is 



e G RB(z = 0) / dE 52 E / dT h G RB(E 5 2,T ; z = 0) = 3.6 x 10 53 ergs Gpc" 3 yr" 1 . (4) 



The average energy release per burst is just the ratio of equations (4) and (3), and is equal to 
^ 2.8 x 10 52 ergs and 8.2 x 10 50 ergs when E$ n = 0.1 and 10" 4 , respectively. This does not 
correspond to the average energy release of detected GRBs, because very energetic bursts are much 
more likely to be detected. One-half of the total energy generated by burst sources comes from 
events with energies = 2.3 x 10 53 ergs. This is a lower limit to the average energy of an event, 
because the use of a single power-law function for Eq in equation (1) does not accurately model 
extremely powerful and very rare events, such as GRB 990123. Consequently, equation (4) is a 
lower limit to the local emissivity determined by fits to BATSE data. This value is not sensitive 
to the choice for E^ m , which is required to be < 0.1 in the statistical study. Thus the local 
volume-averaged GRB energy emissivity is better known than the local GRB event rate density. 

Possible collimation of GRB sources does not alter the energetics arguments made here as it 
does for the event rate calculation, because a smaller beaming fraction is offset by a larger number of 
sources. Neither would beaming affect the efficiency calculations performed below. If GRBs exhibit 
a constant energy reservoir (Frail et al. 2001; Panaitescu and Kumar 2001), the implications for 
GRB source emissivity in the statistical study of Bottcher and Dermer (2000) therefore remains 
unchanged, and in fact implies a distribution of jet opening angles 8 for the sources of GRBs 
(dNcRB/dcosO oc (1 — cos#)~ 1//2 , so that dNcRB/dO oc const when 8 <C 1). Because the rate 
density of GRB sources is inversely proportional to the beaming factor, the radiative signatures 
from a single GRB would be changed due to beaming. Throughout this paper, we quote energy 
emissivities and event rates for uncollimated GRB sources, but additionally consider the constant 
energy reservoir result when calculating the rate of GRBs in the Milky Way. 
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The statistical study of Bottcher and Dernier (2000) is seen to be consistent with other recent 
GRB statistical studies once one recognizes that inefficiencies for generating radiation from the 
GRB event and for detecting emission in the BATSE range have been explicitly taken into account 
in this approach (this point was not considered by Stecker (2000)). Moreover, most burst events 
with T < 100 will not trigger a GRB detector such as BATSE due to the triggering criteria and 
design of burst detectors that have been flown to date (Dermer et al. 1999). These undetected dirty 
fireballs may contribute as much as 50-70% of the total emissivity. 

For example, Schmidt (1999) derives a local emissivity of GRBs in the 10-1000 keV band of 
1.0 x 10 52 ergs Gpc" 3 yr" 1 , which is a factor 36 smaller than the value obtained here. The efficiency 
for the external shock model to produce radiation in the 10-1000 keV band is ~ 5-15%, and ~ 50% 
of the total energy is released in the form of dirty fireballs with Tq < 100 that would not trigger 
BATSE. (The clean fireballs with Tq S> 300 cannot be very numerous.) Insofar as inefficiencies 
for generating 7-ray emission in a colliding shell model are typically 1% or less (Kumar 1999; 
Panaitescu et al. 1999; however, see Beloborodov 2000; Fenimore and Ramirez- Ruiz 1999), and 
that the collision of a relativistic shell with matter at rest allows the greatest fraction of directed 
kinetic energy to be dissipated within the blast wave shell (Piran 1999), we think that equation (4) 
therefore provides a conservative lower estimate for the emissivity of progenitor sources of GRBs 
in the local universe. 

To obtain the emissivity of GRB sources into the Milky Way galaxy, we proceed in two ways. 
The first, following Wijers et al. (1998), is to employ the Schechter luminosity function <&(L)dL 
= ($*/L*)(L/L*) Q ex.p(—L/L*)dL, giving the number density of galaxies with luminosities in the 
range L to L + dL. Assuming that the burst emissivity per galaxy is proportional to the luminosity 
of the galaxy, then 6grb = k J °° dL ■ L ■ &(L), so that k = €GRB[^*L*T(a + 2)] -1 , where T(v) is 
the Gamma function. The energy released by GRB progenitors in a galaxy with luminosity L is 
dE(L)/dt (dE/dVdtdL)/(dN/dVdL) = kL$(L)/$(L), so that 

= g; \ GR *' L v = 2-5 x 10 39 (^) ergs s~\ (5) 
dt <£>*L*T(a + 2) L* v ' 

In the last term of equation (5), we used the results of equation (4) with <!>* = 1.6 x 10 _2 /i 3 Mpc" 3 , 
a = —1.07, and h = 0.65 (Loveday et al. 1992). If the Milky Way is an L* galaxy, then the power 
of GRB sources into the Milky Way is therefore dE/dt > 2.5 x 10 39 ergs s" 
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Scalo and Wheeler (2002) argue that a better approach is to weight the burst emissivity by 
the ratio of the blue luminosity surface density Sl of the Milky Way to the volume-averaged blue 
luminosity density J ga i,B of galaxies in the local universe. Using the expressions Y±l = 20L Q pc~ 2 
and J galyB ^ 1.1 x 1O 8 /iL Mpc" 3 quoted by Scalo and Wheeler (2002), we find dE/dAdt = l.Ox 10 38 
ergs pc -2 yr -1 for the burst power per unit area in the solar neighborhood. For a 15 kpc radius, 
wc then obtain dE/dt = 2 x 10 39 ergs s for the GRB source power in the Milky Way, which is in 
good agreement with the value obtained through the first approach. An advantage of this method 
is to highlight the potentially large variations in the emissivity of GRB sources in different regions 
of a galaxy. 
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The power required to supply the galactic cosmic radiation, assuming that cosmic rays are 
uniformly distributed throughout the disk of the Galaxy, is ~ 5 x 10 40 ergs s _1 (Gaisser 1990). We 
therefore see that GRB sources and the dirty and clean fireballs, collectively referred to as fireball 
transients (FTs), supply a power to the Milky Way that is > 5% of the cosmic ray power, and may 
therefore make an appreciable contribution to cosmic ray production in the Galaxy. The relative 
FT/cosmic-ray power could be much larger if the contribution of clean and dirty fireballs that are 
invisible to GRB detectors (Dermer et al. 1999) is much larger than derived on the basis of the 
single power-law representation of the To- and E'o-distributions. This fraction would also be larger 
if the efficiency for the sources of GRBs to generate 7 rays is smaller than calculated in the external 
shock model used by Bottcher and Dermer (2000). 



2.2. Ultrahigh Energy Cosmic Rays 

The energy density of UHECRs follows from the intensity E 3 dJ/dE = 3.5 x 10 24 eV 2 m~ 2 
gr -i g -i measurec i w ith the Akeno Giant Air Shower Array (Takeda et al. 1998). This expression 
is valid within the experimental error for all cosmic rays with energy 1.2 x 10 19 < £(eV)< 3 x 10 20 
eV, except for being 1.5<r away from the E = 1.4 x 10 20 eV data point. It is accurate to within 
2cr of all data points at 3 x 10 18 < E(eV) < 3 x 10 20 . Above 3 x 10 20 eV, small-number statistics 
dominate. From this expression it follows that the energy density of UHECRs with energy between 
E(eV) and 3 x 10 20 eV is 

uME) (ergs cm-) * ™*^) (* ~ 3 x ^20 eV } ' (6) 

The evidence for a high energy tail above E « 3 x 10 20 eV is unclear due to the small-number 
statistics. Equation (6) should be considered an upper limit to the UHECR energy density, in 
view of the smaller flux of > 10 20 eV particles measured with the monocular HiRes fluorescence air 
shower experiment (Sommers 2002). 

Ultra-high energy particles lose energy by adiabatic losses in the expanding universe, and by 
photo-hadron and photo-pair production on the cosmic microwave background. The mean energy 
loss length x\ oss (E) due to these processes has been recently recalculated by Stanev et al. (2000). 
The loss length for 10 20 eV protons is about 140 Mpc, and this length is also consistent with their 
calculations of horizon distance within which 50% of the protons survive. The values of xi oss at 
E > 6 x 10 19 eV defines the GZK radius insofar as the energy losses are dominated by photo- 
hadronic processes at these energies. The quantity x\ oss (E)/c defines a characteristic survival time 
for particles with energy E. The volume-averaged rate at which astronomical sources produce 
> 10 20 eV particles in the local universe is therefore < 2.4 x 10 -21 ergs cm~ 3 /(140 Mpc/c) = 
1.5 x 10 53 ergs Gpc~ 3 yr -1 , provided that UHECRs traverse roughly straight-line trajectories 
through intergalactic space. If UHECRs diffuse in the intergalactic magnetic field, then the source 
volume contributing to locally observed cosmic rays could be reduced, though with no significant 
change in the local UHECR intensity due to the better trapping of UHECRs in our vicinity. 
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This value is > 2.5 times smaller than the emissivity given in equation (4), so that in principle 
there is a sufficient amount of energy available in the sources of GRBs to power the UHECRs 
(Vietri 1995; Waxman 1995). The conversion of the initial energy of a fireball into UHECRs must, 
however, be very efficient. If nonthermal power-law distributions of particles are accelerated in 
the blast wave, as expected in simple treatments of Fermi acceleration, then hard spectra with a 
nonthermal particle injection index p < 2 place a large fraction of the nonthermal energy in the 
form of the highest energy particles. A large fraction of the blast-wave energy can be dissipated as 
UHECRs even if p > 2 if particle acceleration is sufficiently rapid that particles reach ultra-high 
energies and diffusively escape on the deceleration time scale (Dermer and Humi 2001). 



3. Photomeson and Neutral Particle Production in GRB Blast Waves 

3.1. Photopion Cross Section and Production Spectra 

Only the photomeson process is considered in detail in this paper; photopair and secondary 
production losses involving nucleon-nucleon collisions can be shown to much less important in 
comparison to photomeson losses for ultra- high energy particles in the blast-wave environment. 
The two dominant channels of photomeson production for proton-photon (p + 7) interactions are 
p+j — > p+7T° and p+7 — > n + ir + , which occur with roughly equal cross sections. In the latter case, 
the neutron decays with a lifetime t n = 10 3 s through the /3-decay reaction n — > p + e~ + v e . The 
decay of the charged pion produces three neutrinos and a positron through the chain ir + — ► [i + + 
followed by the decay fi + — > e + + u e + i>„. Neutrino production from photomeson interactions in 
GRB blast waves has been considered earlier (Waxman and Bahcall 1997; Vietri 1998a, b; Rachen 
and Mesazaros 1998; Halzen and Hooper 1999; Waxman and Bahcall 2000; Dai and Lu 2001), but 
usually in the context of an internal shock model. Rachen and Mesazaros (1998a) also consider 
neutron production in the internal shock model. 

To treat neutral particle production, we follow the approach of Stecker (1979) (see also Bottcher 
and Dermer (1998)). The cross section is treated in the (5-function approximation. Thus an in- 
teraction takes place if the photon energy in the proton's rest frame equals "f' p e'(l — //') = €a — 
0.35m p /m e = 640, where primes denote quantities in the comoving frame, 7^ is the proton Lorentz 
factor, e represents photon energies in units of the electron rest mass energy, 6a is the energy of 
the A resonance, and n' is the cosine of the angle between the photon and proton directions. The 
differential cross section for the photomeson production of neutrons and neutrinos produced with 
energy E' is approximated as 

= Oo%' - (1 - ^)}S(E' - m a ' p ) • (7) 

The multiplicity Q = 1/2 for neutrons and Q = 3/2 for neutrinos, noting that we are only consider- 
ing the neutrinos formed from ir + decay (and not from the neutron) . The photomeson cross section 
do = 2 x 10~ 28 cm 2 . Each neutrino carries away about 5% of the proton's initial energy, with the 
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7r + -decay positron receiving another 5%. Thus we let rrii = 0.8m p for neutrons and mi = 0.05m p 
for neutrinos, with the units of the proton rest mass m p defining the units of E' . 

The neutral particle production spectrum in the comoving frame is therefore 

N'(E') - I d 7 ' p N^' p ) f~ de>n> ph (e>) 1^ dul{\ - „t) (^gp) , (8) 

where N' p (p( p ) gives the nonthermal proton spectrum in the comoving frame, and n' ph (e')de' is the 
differential number density of soft photons, assumed to be isotropically distributed in the blast-wave 
fluid frame, with photon energies between e' and e' + de'. Substituting equation (7) into equation 
(8) gives 

N[{E') = W r de' e'- 1 n' ph (e') . (9) 

Ztj J mi t A /2E' 

The production spectrum of neutral particles as measured by an observer can be approximately 
obtained by noting that the differential time element dt = dt'/T, and particle energy E = TE' } 
where unprimed quantities refer to observed quantities. Redshift effects are not considered in this 
section. In more accurate treatments, a full angular integration over the production spectrum 
should be performed, which is especially important if the outflow is collimated. In the present 
treatment, it is adequate to use the simpler relations for dt and E. Thus N[{E') = Ni(E), and we 
have 

Ni(E) = N' p {E/Tm t ) I(y) , (10) 

where 

I(y) = 1^ de' e'- 1 n' ph (e') , and (11) 



3.2. Blast Wave Dynamics 

We consider the case of an adiabatic blast wave decelerating in a uniform surrounding medium 
with density no- 1 When r> 1, the blast wave evolves according to the relation 

T(x) = - r ° (12) 

(Chiang and Dermer 1999), where x is the distance of the blast wave from the explosion center, 
and the deceleration radius 

Xd ee ( r f° 2 )V3 ^ 2.1 x 10 16 (-#^)V3 cm (13) 



1 This treatment is reasonably consistent with the statistical treatment of GRBs by Bottcher and Dermer (2000) 
There the radiative regime that provided the best fit to the GRB statistics is nearly adiabatic, with the blast wave 
decelerating as f a x~ 1,7 , compared to F oc a;" 1,5 in the fully adiabatic limit. 
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(Meszaros and Rees 1993), where I^oo = To/300. The rate at which nonthermal proton kinetic 
energy is swept-up in the comoving frame of an uncollimated blast wave is 

E' ke = 4TTX 2 n f3c{m p c 2 )r(T - 1) (14) 

(Blandford and McKee 1976), where (3 = y/l — T~ 2 . Thus the accumulated nonthermal kinetic 
energy at radius x is 

E> (x) - ( X dx\ dt '\E' - E ° J ' X<<Xd , (15) 

where the largest value of x in the second asymptote stems from the T S> 1 restriction, and 
dx = f3Tcdt'. 

It is convenient to relate the observer's time t to x, and describe blast-wave evolution in terms 
of the dimensionless time r = t/t^ where the deceleration timescale (for the observer) is 

t d = (1 + z)$- = 7.7(1 + z)(-^) 1 / 3 s (16) 

(Rees and Meszaros 1992; Meszaros and Rees 1993). Because dt = dx/T 2 c, 



Xd 

Likewise, 



r , for r <C 1 r 8/3 

(4r)V4, forl<<r<<Tf " TTFVW' fOTr<<ro " 



1, forr<l ^ 1 

= \ (4r)~ 3 / 8 , forl«T«Tf " Vl + (4r)3/4 ' 



E,\\t\ forr«l ^ E Q ,^ 



and 

eL{t > a r, { (4,,3/ 8 , i ; ; « & - ^ra^ • f « - « ^ ■ < i9 > 

The expressions on the right-hand-sides of equations (17)-(19) accurately bridge the early and late 
time behaviors of the asymptotes. 



3.3. Comoving Proton, Electron, and Photon Spectra 

A power-law distribution of nonthermal protons with number index p is assumed to be accel- 
erated in the blast-wave. Because protons and ions are swept up with Lorentz factor T and are 
then subsequently accelerated, we represent the nonthermal proton distribution by the expression 

T ) = KrT p ELi t P , l'v V > f - r(r) < 7 ; < 7 max (20) 
m p c 2 (T 2 -P - 7ma/) 
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(Bottcher and Dermer 1998). The term £ represents the fraction of swept-up particle kinetic energy 
that is transformed into the energy in the nonthermal proton distribution and could, in principle, 
be as large as ~ 0.5. Not more than ~ 10-20% of the total nonthermal proton energy could, 
however, be radiated if the treatment is to remain consistent with the assumption of an adiabatic 
blast wave. The term 7 r ' nax , giving the maximum proton Lorentz factor in the blast wave frame, 
must be > 10 10 for GRBs to account for UHECRs. Rachen and Mesazaros (1998) define limits on 
various acceleration scenarios that give large values of 7max- Particle spectra from gyroresonant 
acceleration due to pitch-angle scatterings and stochastic energy diffusion in particle interactions 
with plasma waves can give 7m ax ~ 10 10 an d can produce nonthermal spectra with p > 1, though 
the exact value depends on the spectrum of the turbulence (Schlickeiser and Dermer 2000; Dermer 
and Humi 2001). 

A nonthermal electron spectrum is also assumed to be accelerated in the blast wave with 
the same index p as the nonthermal protons. Following Sari et al. (1998) (see also Dermer et al. 
(2000a)), we represent the nonthermal electron spectrum by the expression 

AM) . (. - W { p-,^ £ I (2D 

where N e = 47rx 3 no/3 is the total number of swept-up nonthermal electrons and is the electron 
Lorentz factor. In the slow cooling limit, 70 = j m , 71 = j c , and the steady-state electron spectral 
index s = p, whereas in the fast cooling limit 70 = j c , 71 = 7 m , and s = 2. Here the minimum 
electron Lorentz factor j m = e e (p — 2)Tm p /[{p — l)m e ] and the cooling electron Lorentz factor 
7 C = 3m e /(16m p e£noc<7TT 3 i), where e e and e# are parameters describing the swept-up kinetic 
energy transferred to the electrons and the magnetic field, respectively (Sari et al. 1998). The 
magnetic field B is defined through the expression 

B = ^32Tm mpC 2 e B r(r - 1) 9* 0.39^/e^ T G . (22) 

We let 72 = 4 x 10 7 e max /[i3(G)] 1 / 2 (Chiang and Dermer 1999) and take e max = 1 in this paper. 

We consider only nonthermal synchrotron emission here. Synchrotron self-absorption and syn- 
chrotron self-Compton (SSC) processes are treated by Dermer et al. (2000a), including a comparison 
of the analytic results to detailed numerical simulations. Given the parameters used here, the ne- 
glect of synchrotron self-absorption is not important for photomeson production, but the inclusion 
of Compton processes could, however, depress the intensity of the low-energy photon spectrum 
when ee -C e e . 

In the 5-function approximation for the synchrotron emissivity, the photon production spec- 
trum 

W = l^'-^e-^Ki^) , (23) 

where v$ = vq{B) = 4ccjt'Ub/3, ub = B 2 / (87rm e c 2 ) is the magnetic-field energy density in units of 
m e c 2 cm -3 , and en = B/B cr = 5/4.413 x 10 13 G. The magnetic field is assumed to be randomly 
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oriented. This formula is accurate to better than a factor-of-2 except near the endpoints of the 
distribution (see Fig. 2 in Dermer et al. (2000a)). By substituting equation (21) into equation (23), 
we obtain the comoving photon density 

ilo 2/3 (e'/e H )-^ 

7<rV/^r (s+1)/2 





n' h (e') = K 



„ s -i^+i- S(e / /eff) -(p+2)/2 



for e'/e H < 7o 
for 7 2 < e'/e H < 7l 2 
for 7l 2 < e'/e H < 7 | 
for 7 | < e'/e H , 



(24) 



where 



8 4irx 2 ce 2 H 



487rm e c 2 



(25) 



Substituting equation (24) into equation (11), performing the integrals, and defining 7flj j = 
Te A /(2e H -f?) for i = 0, 1, and 2, we obtain 



fo, 



7(7) = 2K < 



_ _ , -(P+2J/2 _ ^-P- 2 1 

-2e w7 V '2 J 

( P +2)- i 7o - i 7 r i - s ( 7 r p - 2 -7 2 -^ 2 ) 

+(,+i)- i 7o s - i [(g^)-( s+i )/ 2 - 7 r- 1 
( P +2)- i 7o - i 7 r i - s ( 7 r p - 2 - 72 - p - 2 



for 7 < 7„ i2 

for 7„ j2 < 7 < 7n,i 

for 7n,i < 7 < 7n,0 



(26) 



+ ( S + 1)- 1 7 S - 1 

-4/3 



] , for 7 > 7„ i0 , 



where 7 = E/rrii. The (reA/2e#7) term dominates each of the branches of equation (26). A good 
approximation to 1(7) is therefore 

, for 7 < 7„ j2 

(p + 2)-H^) s+1 t^)^ 2 , for 7 n,2 < 7 < 7n,l 



2K 



(a+ i)-i (^)^)/ 2 , 

3 / 7 \2/3 
l> 4 V7„, ^ ' 



for 7„,i < 7 < 7„,o 
for 7 > 7n,o • 



(27) 



The production spectrum Ni(E) of neutral particles formed through photomeson production is 
therefore given by equation (10), but with I(y) replaced by either 1(7) or I a p(l) given by equations 
(26) or (27), respectively. 



3.4. Energy-Loss Timescales for High Energy Protons 

Energy-loss timescales are derived in the comoving frame for protons that would have energies 
E as measured in the observer frame. These timescales are compared with the comoving time t' 
passing since the initial explosion event; clearly if the energy-loss timescale is long compared with 
the available comoving time, then only a small fraction of the particle energy can be extracted 
through that process. From the relation dx = flTcdt', we obtain the comoving time 

4i ~ r + / T ' for r < 1 T t d r 8/3 . . 

* = { 2 (4r)5 /8 > for 1<<T<< r 8/3 ^ , for r « Iq . (28) 
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The photopion energy-loss rate is 

1 c f°° f 1 
t'pj - 5 2 _/ de ' n ph( e ') J i ~ /*') cyy-^n^V) ' ( 29 ) 

(compare eq. [8]), where the factor 1/5 takes into account that « 5 interactions are required 
for a high-energy proton to lose a significant amount of its energy. Here we consider both the 
P7 — ► 7r + n and p7 — ► 7r°p chains, because both will compete against other energy-loss processes. 
Thus a pl ^ n (e',n') = a 5[fi f - (1 - e A /7 p e')}, giving 



where 7 = E/m p = Tj' p . Hence 



# a ^ /<7) ' <30) 



-(#)„ - ^ ^) • (3D 



The energy- loss rate through photopair (p + 7 — ► p + e + + e~) production is small compared to the 
photomeson energy-loss rate at very high energies because of the greater energy loss per scattering 
event in photomeson production. Although photopair production could dominate the energy-loss 
rate for protons with 7^ <C 10 8 , it is not important for the highest energy protons and is not treated 
here. 

The proton synchrotron loss rate is given by 

g v oi P _ eBno(cm- 3 )r 2 7 ; 
p ' syn {m p /m e f 3.2 x 10 19 s ' 1 j 

The importance of this process for producing high-energy 7 rays from GRB blast waves has been 
considered by Vietri (1997) and Bottcher and Dermer (1998). 

The secondary production rate is t 1 ^ = n'a pp c, where n' = 4Tuq from the shock jump condi- 
tions, and a pp = 30 mb. The secondary production efficiency r] pp = t/t pp = i7 pp t/[1-|-(4t) 3 / 4 ], where 
f] pp = &noo pp Xd = 5 x 10 _9 £ , 52 3 (no/r3oo) 2/ ' 3 . The secondary production efficiency increases through 
the afterglow phase, and reaches a value of rj pp = 1CT 7 'E^n^ 3 at r = T 8 / 3 . Unless no » 10 8 cm -3 , 
the efficiency for this process will be low (see Pohl and Schlickeiser (2000) and Schuster, Pohl, & 
Schlickeiser (2002) for treatments of this process in GRBs and blazars, respectively). 

Fig. 1 shows results of calculations of the ratio 77 of the comoving time to timescales for pho- 
tomeson production (open circles) and proton synchrotron radiation (filled circles). The timescales 
are calculated at different observer times as a function of observed proton energy E, up to the 
maximum proton energy defined by the synchrotron radiation limit given in Appendix A. The 
chosen parameters in Figs, la and lb are typical of those used to fit GRBs in the prompt and 
afterglow phase, respectively, and are listed in Table 1. In both cases, we use a total energy release 
Eq = 2 x 10 53 ergs, which is near the mean value of the energy release distribution (see Section 
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2.1). The value p = 2.2 is similar to that deduced in fits to afterglow GRB spectra of GRB 990510 
(Harrison et al. 1999) and GRB 970508 (Wijers and Galama 1999); a value of p much steeper than 
~ 2.2 will make the energetics of UHECR production problematic. 

Other than Eq and p, Fig. la employs the parameter set in Fig. 1 of Dermer et al. (2000b) which 
was shown to give good fits to burst spectra during the 7-ray luminous phase of GRBs (Chiang 
and Dermer (1999); there we used Eq = 10 54 ergs and p = 2.5). The remaining parameters used 
in Fig. la are Tq = 300, tiq = 100 cm -3 , es = 10 -4 , and e e = 0.5. Even with such a large value 
of e e , the blast wave evolves in the adiabatic limit because the electrons are in the weakly cooling 
regime. We also take £ = 0.5. The dotted lines show the photopion timescales obtained using 
the approximate expression for I a p{l) in equation (27). Fig. lb uses parameters that are typical 
of those used to model the afterglow spectra of GRBs (Harrison et al. 1999; Wijers and Galama 
1999), and are the same as Fig. la except that cb = 0.1 and e e = 0.1. The latter choice ensures 
that the GRB blast- wave evolution is nearly (though not quite; see Bottcher and Dermer (2000a)) 
adiabatic, given that the electrons are strongly cooled during the prompt phase and much of the 
afterglow phase for this larger value of e#. The major difference between the fits derived to GRB 
spectra during the prompt and afterglow phases is thus the stronger value of field at later times. 
Other arguments that the magnetic field evolves to its equipartition value following the prompt 
phase are given by Dermer et al. (2000b). 

As can be seen from Fig. 1, the relative timescales for photomeson production in the external 
shock model usually dominates the other processes, and approaches or exceeds unity for the highest 
energy protons during the afterglow phase. Thus a large fraction of the energy contained in the 
highest energy protons is converted into an internal electromagnetic cascade and lost as photomeson 
neutral secondaries. The largest proton energies are constrained by the synchrotron limit given in 
Appendix A, but still exceeds 10 20 eV, in accord with the hypothesis that UHECRs are accelerated 
by GRBs. Protons with observed energies > 10 18 eV therefore lose a significant fraction of their 
kinetic energy through photomeson production which is transformed into neutrons, neutrinos, and 
high energy leptons. The leptons generate high energy gamma rays during an electromagnetic 
cascade in the blast wave (Bottcher and Dermer 1998). When e# > 0.1, proton synchrotron 
losses can dominate photomeson losses during certain phases of the evolution. Although secondary 
production can be the dominant proton energy loss process at E < 10 16 eV, its importance is 
negligible unless the CBM density is very high. 

When the relative timescales exceed unity, a large fraction of the proton energy is radiated 
away during the comoving time t', and the proton distribution will strongly evolve through radia- 
tive cooling. When this occurs, a thick-target calculation is required to calculate total neutrino 
and neutron emissivity. This regime begins to be encountered here, but a complete treatment of 
photopion production will require solving a transport equation that is beyond the scope of this 
paper. We also note that the efficiency ratio rj = t' /t' (process) is actually shorter when the injec- 
tion index p becomes larger, because the energy density of the soft photons is then concentrated 
into a narrower bandwidth and is therefore more intense. Nevertheless, much less energy of the 
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total GRB energy is radiated through photomeson production when p » 2, because the total GRB 
energy carried by the highest energy protons is much smaller. 



3.5. Instantaneous and Time-Integrated Production Spectra 

It is simple to derive the characteristic spectral behavior of neutrons or neutrinos produced 
in the external shock model. Using equations (27) and (20) in equation (10), we find that the 
instantaneous production spectra, multiplied by E 2 , follow the behavior 

'0, for E/m,i < max(r,7„ i2 ) 

E (4- P )/2 j for max(rj ln 2) < E/m . < 7nji 

E 2 Ni{E) oc { E$+ s - 2 p)I 2 , for 7n>1 < E/m t < 7n , (33) 

E -p+5/3 j for 7n () < £ /m . < 7L ^ 

, for E/m; > 7^ 



j,max 

,,max i 



where ^L,max is given by the synchrotron radiation limit, equation (Al). These spectral indices 
are two units larger than particle injection number indices. It is also assumed in these expressions 
that T < 7 n> i and "fL,max > 7n,o> but it is simple to generalize the results when this is not the case. 
The instantaneous production spectra are very hard at low energies, with Ni(E) oc E~' p l 2 £S E^ 1 
when p ~ 2. Irrespective of whether we are in the fast cooling (s = 2) or slow cooling (s = p) 
regime, the spectra soften to Ni(E) £S E~^/ 2 for p ~ 2, although the spectra still rise in an E 2 Ni(E) 
representation. At energies E > ruijug, Ni(E) oc E 1 ^ -1 / 3 £ E~ 7 / 3 , where the —7/3 behavior holds 
when p ~ 2. The E 2 Ni(E) peak energy is carried primarily by particles with energy 

^ _ mJtA _ mJ 2 6 A (r/300) 2 J 2 x 10 17 eV , for neutrons , . 

pk ~ man '° ~ 2e Hl 2 ~ 2e br ~ (e br /0.1) \ 10 16 eV , for neutrinos . { ' 

In this expression, the break energy e br is the photon energy separating the e~ 2 / 3 portion of the 
synchrotron emissivity spectrum produced by an electron distribution with a low-energy cutoff 
from the higher-energy portion of the synchrotron spectrum. As is well known, this often occurs at 
energies ~ 50 keV - several MeV during the prompt phase of GRBs (Cohen et al. 1997). Equation 
(34) also follows from elementary considerations. 

Figs. 2a and 2b show instantaneous production spectra at different observing times for neutrons 
and neutrinos, respectively. Fig. 2a employs the parameter set (A) for the prompt phase of GRBs 
and Fig. 2b uses the parameter set (B) that better represents afterglow data (see Table 1). The 
peak of the E 2 Ni(E) spectrum at Ep^, above which Ni(E) oc E~ 2 ' 53 , occurs at early times in the 
instantaneous spectra of Fig. 2a. The transition to the I? -p ~ 1 / 3 portion of the spectrum is not seen 
after ~ 10 s is Fig. 2a nor in the Fig. 2b spectra. This is because e^r reaches such low energies that 
E p k would occur above the maximum energy defined by equation (Al). The maximum energies of 
the neutrons reach or exceed ~ 10 21 eV, but the maximum neutrino energies only reach < 10 20 eV 
due to the smaller amount of energy transferred to each neutrino in the photomeson production 
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process (compare cross section (7)). The production spectrum breaks from Ni{E) oc E~ - 9 at low 
energies to Ni(E) oc E~ 1 ^ at intermediate energies in Fig. 2a, because the electrons distribution 
starts to evolve in the uncooled regime. In contrast, the spectrum above the break in Fig. 2b is 
slightly softer with Ni(E) oc E~ ir ' ', because the electron distribution evolves in the strongly cooled 
regime. 

We also show the time-integrated production spectra of both the neutrons and neutrinos for 
parameter sets (A) and (B) in Figs. 2a and 2b, respectively. Here we integrate the instantaneous 
production spectra over all times until the blast wave reaches x = x^T^ 3 , where it has decelerated 
to mildly relativistic speeds. The time-integrated spectra retains its Ni{E) oc — p/2 behavior at 
10 12 < E(eV) <C 10 17 eV. For the prompt-phase parameter set (A), the time-integrated spectra 
steepens to a Ni(E) ,25 E -1 ' 8 behavior above the value of E p k evaluated at t = tj, and then cuts 
off at a maximum energy determined by equation (Al) at r = 1. For parameter set (B), the time- 
integrated spectrum remains very hard, with Ni(E) £5 E~ p / 2 , up to nearly the maximum energy 
defined by ma L>max . 

The time-integrated spectra in Fig. 2 imply both the total energy release and the energies of 
the produced neutrons and neutrinos that carry the bulk of this energy. Neutrons with energies 
between ~ 10 18 and ~ 10 21 eV carry ~ 10 51 ergs of energy for the chosen parameters. Neutrinos 
carry ~ 3/20 as much total energy as the neutrons in an energy range that is ~ 20 times smaller 
than that of the neutrons. The ratio of the energy carried by either neutrons or neutrinos to the 
total explosion energy Eq, here called the production efficiency, is therefore ~ 1% for neutrons and 
~ 0.1% for neutrinos. Fig. 3 shows calculations for the neutron and neutrino production efficiencies 
as a function of £54. The neutron production efficiency increases with increasing Eq and reaches 
a few per cent when E54 = 1. Parameter set (A) gives better efficiency at large values of Eq 
than set (B), but poorer efficiencies when £54 < 0.2. The production efficiency is only weakly 
dependent upon Tq, but depends strongly upon p, as outlined earlier and shown in the inset. The 
maximum efficiency occurs when p ~ 2.1. According to the statistical treatment of the external 
shock model described in Section 2.1, ~ 50% of the total GRB energy is radiated by explosions 
with Eq > 2 x 10 53 ergs. Thus we find that > 1% and > 0.2% of this energy is converted into 
high-energy neutrons and neutrinos, respectively, if the UHECR/GRB hypothesis is correct. This 
will have the observable consequences described in Sections 4 and 5. 

The efficiencies for neutral particle production correspond to GRBs with smooth profiles. In 
the external shock model result, smooth-profile GRBs result from blast wave deceleration in a 
uniform surrounding medium. We have chosen uq = 100 cm -3 . Because rj oc y/fiQ (Appendix B), 
smooth profile GRBs could produce neutrinos and neutrons with even smaller efficiencies if no < 1 
cm -3 . Interactions with an inhomogeneous and clumpy CBM are thought to produce the short 
timescale variability observed in rapidly variable GRBs in the external shock model (Dermer and 
Mitman 1999; Dermer and Bottcher 2002). Under these circumstances, neutrino production could 
be considerably enhanced. Thus the production calculations only apply to GRBs which display 
smooth profiles. 
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3.6. Temporal Behavior of Production Spectra 

The temporal indices of the particles formed through photomeson production can be obtained 
by examining equations (10), (20) and (27), noting the temporal dependences of the various terms. 
Writing equation (10) in more detail, we have 

A>/^_ C^oeAr (p-2)£E% e (r) E 

N ^ E ) ~ Kp ^7 ( TW [Iap(E/mi)/K\ . (35) 

m p C 2 (T 2 p -7max ) im « 

provided T < E/rrii < ^y' max . The coefficient K has been extracted from the I a p(l) term, and varies 
according to K(t) oc x (equation (25)), so that it has the time dependence given by equation (17). 
The time dependences of V and E' ke are given by equations (18) and (19), respectively. It then 
becomes necessary to determine the time dependences of = Te/\/2eH^Q and therefore of the 
that enter into equation (27), noting that eh oc B oc V. 

The temporal behavior of the neutron and neutrino production time profiles, or "light curves," 
depends on whether the electrons are in the slow or fast cooling regimes (Sari et al. 1998). Because 
of the progressive weakening of the magnetic field in the standard blast- wave model, the fast cooling 
regime will exist only if the cooling electron Lorentz factor 7 C is less than the minimum electron 
injection Lorentz factor j m at r 1. Using the expressions for j c and j m following equation (21), 
we therefore find that the nonthermal electrons will evolve in the fast cooling regime at least during 
some stage of the blast-wave evolution if 



r n > f r 



0.16 . 1 ,p 



1 / 2 E V 1 A e eeB p-2 



(^4)] 3/4 , (36) 



using equation (16) for t^. When equation (36) does not hold, the system is always in the slow 
cooling regime. For example, if we vary only Tq in parameter set (A), there will be some evolution 
in the fast cooling regime when Tq > 70. There will be evolution in the fast cooling regime for 
essentially all values of To 3> 1 with parameter set (B). Note that the baryon- loading factor Tq 

- — 1/2 

separating the different cooling regimes is quite sensitive to no, with T oc n Q 

Fig. 4 is a sketch of the temporal indices of particles produced with different energies as a 
function of dimensionless time r. First consider the outside boundaries of the temporal-index 
plane. Particles will only be produced if E/rrii > T. This defines the lower region bordered by 
the short dashed lines. Due to threshold effects, another limit to low-energy particle production 
arises from photomeson threshold effects due to the upper cutoff of the highest energy photons at 
e' ~ ^2 e H- Neutrons and neutrinos will not be produced with energies mj7 ni 2 < m^t^l{2^\tu) 
because of this cutoff. For the synchrotron radiation limit given by equation (Al) — but now 
for electrons — 72 oc B 1 / 2 , so 7^2 oc as sketched by the dashed-triple-dotted lines. The 
synchrotron radiation limit for protons defines the upper boundary for the highest energy particles 
that are produced. It is oc T 1 / 2 , and is shown by the dot-dashed lines. Different Fermi acceleration 
models could give different maximum particle energies, but all would likely be bounded from above 
by this limit due to the competition between the synchrotron loss and acceleration rates. 
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The blast-wave system, as illustrated in Fig. 4, passes through a fast cooling regime. Thus there 
are two dimensionless times r< and r> defined by the relation j c = 7 m that bound the period when 
the blast wave is in this regime. The blast wave is in the slow cooling regime when either equation 
(36) fails to hold or, if not, when r < r< and r > r>. We define 7„ im = r^A/ '{^Hlm) ^ ^ 2 
and 7„ iC = reA/(2e#7 2 ) « T 6 t 2 . Hence 7 n , m oc 1 and 7„ iC oc r 2 for r < 1, and 7„ im oc r 3/4 and 
7n,c oc r~ 1//4 when 1 < r < r 8 / 3 . The behaviors of j n ,m and 7„ iC are indicated by the thick lines 
and the double lines, respectively, in Fig. 4. 

In the slow-cooling regime, 70 = 7 m and 71 = j c . In the fast-cooling regime, 70 = 7 C and 
7i = 7m- It is straightforward though tedious to derive the temporal indices x displayed in Fig. 4 
using the above relations and equation (27) in equation (35). The important point to notice is how 
hard the values of x are m the afterglow phase. The highest energy particles with 7 ni o < 7 < r7^ ax 
are due to interactions with the e~ 2 / 3 part of the soft photon spectrum. In the uncooled regime, 
X ~ —0.25 when p ~ 2, so that the bulk of the energy is radiated at late times. Because of the rapid 
decay of with time, however, this phase does not persist very long. Nevertheless, the temporal 
indices in the lower and intermediate energy regimes are % = —0.75 and x — —0.875, respectively, in 
the afterglow phase when p ~ 2. Thus the bulk of the energy is still radiated at late times. Although 
this temporal behavior will be difficult to detect from neutrinos and neutrons from GRBs, they are 
relevant to the high-energy gamma-ray spectrum observed from GRBs. Charged pions will decay 
into leptons, which will scatter soft photons to high energies to generate a cascade, and neutral 
pions from p + 7 — ► p + 7r° will decay to form 7 rays that can pair produce until the photons are at 
sufficiently low energies to escape. As noted by Bottcher and Dermer (1998), the temporal decay 
of the high-energy emission from hadrons is much slower than the synchrotron decay. Thus high- 
quality GeV observations of GRBs could reveal the presence of a high-energy hadronic component, 
though it must be carefully distinguished from the SSC component, for example, by its spectral 
characteristics. The Gamma ray Large Area Space Telescope (GLAST) mission 2 will be well-suited 
to measure the 7-ray afterglow of GRBs and thus test for an energetic hadronic component in GRB 
blast waves. 

Fig. 5 shows calculations of the neutron and neutrino production time profiles at 10 12 , 10 15 , 
10 18 , and 10 20 eV for parameter sets (A) and (B), as described in the figure caption. Here we have 
multiplied the E 2 N(E) spectra by observing time t in order to reveal the time during which the 
bulk of the energy is radiated. For these parameters, roughly equal energy is radiated per decade 
of time, except at the very highest energies and during early times. From the analytic results for 
the temporal index in the afterglow phase, we find that x — —0.32 and % = — 1.1 at intermediate 
energies, in agreement with the calculations. The lower energy regime with x — —0.975 is not 
encountered here. The abrupt cutoffs at early and late times are due to the definite ranges of 
particle energies implied by the analysis. 



2 http://glast. gsfc.nasa.gov 



-20- 



4. Neutrinos from GRBs 

The detailed calculations provide a lower limit to the neutrino fluxes if the UHECR/GRB 
hypothesis is correct. Even within the context of the external shock model, other effects could 
enhance the neutrino emissivity. For example, reverse shock emission provides additional soft 
photons that would enhance photomeson production (Dai and Lu 2001). Larger neutrino fluxes 
could also be obtained if we relax the assumption that the surrounding medium is uniform, which 
is probably the case in many GRBs, in view of the short timescale variability observed in their 
7-ray emission (Dermer and Mitman 1999). Also important is the uncertainty in determining the 
rate density of dirty and clean fireballs. 

Before displaying calculations, it is useful to make an estimate of the neutrino background 
expected from GRBs. The energy density of high-energy neutrinos from GRBs is 

u u ~ r/ u - E(l) • e GRB (z = 0)-f H ~6x 10 _1 V ergs cm" 3 (37) 

where we use equations (4) and (2) to give the mean emissivity at z ~ 1, and let the Hubble 
time tu = 10 10 yr. The term rj u represents the production efficiency which, as we have seen, is 
~ 0.1%. For p ~ 2, the principal behavior of the time-integrated GRB neutrino spectrum varies 
oc E^ 2 ~ E~ x up to some maximum energy -E^ max ~ 10 18 -10 19 eV (see Fig. 2 and App. B), so that 
the diffuse neutrino background flux <& V {E V ) also varies oc E~ x . Using equation (37) to normalize 
this flux, we find 

E ■ * (E ) ~ 9xlQ "' 9 ^ ir3 ) cm" 2 s- 1 sr- 1 (38) 
U A (iW/10 18 eV) Cm S 81 • ^ 

GRBs will therefore produce a diffuse neutrino flux at the level of 10~ 9 GeV cm" 2 s" 1 sr" 1 at 
E v > 10 18 eV. This flux 3 is comparable to other estimates of cosmological neutrinos above ~ 10 16 
eV (Gaisser et al. 1995; Yoshida and Teshima 1993; see also Stecker et al. 1991). 



4.1. Diffuse Neutrino Background 

Each GRB produces a time- integrated neutrino spectrum (dE/dEf, m ) =Ef, m (dN/dE° m ), where 
is the energy of the emitted neutrino. The luminosity distance di is defined so that the 
relationship dE/dAdt = {AiTd 2 L )~ l {dE / 'dt em ) holds, where dt em is the differential element of time 
in the emitter frame. Because dt = (1 + z)dt om and E v = Ef, m /(l + z), we find that 

dE _ (1 + z) 2 dE 
dAdE v Audi dE™>' { ' 



3 A similar estimate for the diffuse extragalactic 7-ray background, though now assuming a spectrum varying cc e^ 2 
as might be expected from an electromagnetic cascade in the GRB blast wave, gives a diffuse flux of e 2 • "J> 7 ~ 0.1?7 7 
keV cm" 2 s" 1 sr" 1 . The 7-ray production efficiency 7/ 7 could reach 10%, but this emission still falls well below the 
observed diffuse extragalactic 7-ray background, which has an intensity of ~ 1 keV cm" 2 s" 1 sr" 1 between ~ 100 
MeV and ~ 100 GeV (Sreekumar et al. 1998). 
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The differential event rate observed from bursting sources with comoving density n is dN = (1 + 
z)~^ncd 2 L d£l x\dt cm /dz\dz (e.g., Weinberg 1972; Dermer 1992). For a Friedmann- Robertson- Walker 
universe, \dt cm /dz\~ l = H (l + z)^(l + Sl m z)(l + z) 2 - n A (2z + z 2 ) (Totani 1999). The diffuse 
flux of neutrinos produced by the superposition of GRBs throughout the universe is therefore 

d E c r dz r dTo r dEo n^^z).^^)^ 



dAdE u dtdn AttHq J J 1 J (1 + z) V(i + M(l + z ) 2 - ^a(2^ + z 2 ) 

Calculations of the diffuse neutrino background, using equation (1) with standard parameter 
sets (A) and (B) in equation (40), are shown in Fig. 6. As before, we use a cosmology with 
(^o,^a) = (0.3,0.7) and h = 0.65. The relevant units conversion is c/(47rF -Gpc 3 -yr) = 1-22 x 
10 -63 cm~ 2 s _1 . The solid curve is for parameter set (A), and the dotted curve is for parameter set 
(B). The calculations are at the level of E v ■ §„{E V ) ~ 10~ 18 cm" 2 s" 1 sr" 1 at E v <C 10 18 eV. Thus 
the estimate of the diffuse neutrino background (38) is in accord with these results. Parameter set 
(A) produces a more luminous neutrino flux at lower energies because there are greater number of 
high energy soft photons, due to the smaller magnetic field used in this parameter set. It is not 
clear in this representation, but there is approximately equal energy fluxes in parameter sets (A) 
and (B), but most of the energy for set (B) is carried by neutrinos with energies between 10 18 and 
10 19 eV. 



4.2. GRB Neutrino Event Rate 

The number of neutrino events that would be detected per year by a muon detector with an 
effective area of A(km 2 ) due to upward-going neutrinos is 

f°° d W dF 

w — *■ - 3M y 107 ' 1010 A 2 '1 ~t' JME^m ■ ■ < 41 > 

In this expression, P V ^^(E V ) is the probability that a neutrino with energy E u , on a trajetory 
passing through a detector, produces a muon above threshold. From the work of Gaisser and Grillo 
(1987) and Lipari and Stanev (1991) as summarized in Gaisser et al. (1995), we use the following 
approximation to calculate neutrino event rates: 

( 5.2 x 10" 33 [^(eV)] 2 - 2 , for 10 9 < ^(eV) < 10 12 
P V ^(E V ) w \ 3.3 x 10- 16 [^(eV)]°- 8 , for 10 12 < £„(eV) < 1.2 x 10 15 (42) 
[ 1.1 x 10- n [^(eV)]°- 5 , for 1.2 x 10 15 < ^(eV) 

(see also Dai and Lu (2001)). 

The number of events detected from a GRB at redshift z is 

(l + z) 2 10 w A /*°° dE v , dE 



v events 



4vr4 



[°°dE dE 

Jo ^W W ^" 



(43) 
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Here, in contrast to equation (40), it is necessary to evaluate d>L explicitly. It is given by 

a±»T , (44) 



The size distribution of neutrino events can be obtained by evaluating the quantity 

N v (> Events) = -5- dz dF dE , (45) 

-no Jo Jl Jo (1 + zY^[l + U m z)(l + z) 2 - U A (2z + z 2 ) 

where a contribution to the integral occurs only if the number of neutrino events, calculated through 
equation (43), exceeds N events . 

This calculation is displayed in Fig. 7. Neutrino detectors at energies S> 10 12 eV are more 
sensitive to neutrino number flux rather than energy flux, so these neutrino spectra do not regret- 
tably yield large numbers of neutrino events per year. Very weak neutrino fluxes are predicted 
in the external shock model of GRBs for smooth profile GRBs, and we predict that no neutrinos 
will be detected in coincidence with such GRBs. Detection by km 2 detectors of multiple neutrino 
events from smooth-profile GRBs would probably rule out the external shock model. Obversely, 
the lack of detection of neutrino events from GRBs with smooth profiles is fully consistent with the 
underlying assumptions of this study. The predicted neutrino signal from highly variable GRBs will 
require more study. Detectable neutrino emission with km 3 neutrino telescopes is predicted within 
an internal shell model (Waxman and Bahcall 1997), though it is not clear if these predictions also 
apply to smooth-profile GRBs. 



5. Radiation Halos from Neutrons Produced by GRBs 

After production, neutrons with Lorentz factor j n = 10 10 7io will travel a characteristic distance 
\ n = cj n t n = IOO710 kpc before their numbers are depleted by /3-decay. Figs. 2a and 2b show that 
for strong GRBs, > 1% of the explosion energy is carried by neutrons with E ~ 10 18 -10 20 eV, or 
0.1 < 710 < 10. Thus a halo of neutron-decay electrons, protons, and neutrinos will be formed 
around the site of a GRB that extends over a size scale of ~ 10 kpc - 1 Mpc. The characteristic 
neutron-decay lifetime is t n ^ n ~ 3 x 10 5 7io yrs. If the /3-decay electrons radiate on a timescale that 
is much shorter than the neutron-decay timescale, then the peak power of a single energetic GRB 
explosion in the extended nonthermal radiation halo could reach 

d-^haio „ m t- m e lO 54 ^ 36 -E54 _i 

— ; — ~ 0.01JF ~ 10 J- ergs s , (46) 

dt m p 7„t„ 710 

where T = 0.1 is a temporal correction factor. This assumes a neutron-production efficiency of 1%, 
as applies to GRBs with Eq > 2 x 10 53 ergs (see Fig. 3). The power in the outflowing neutron- 
decay protons is ~ m p /m e larger than that given by equation (46), but will only be detected if 
the protons can radiate this energy or deposit it in a galaxy halo to be radiated through other 



-23- 



processes. Given the rarity of powerful GRB events, it is unlikely that the nonthermal halo from 
neutron-decay electrons surrounding a galaxy will consist of the superposition from several GRBs 
at the frequency of greatest luminosity, but will instead be formed by a single event. If GRBs are 
stongly beamed, however, the radiation halos will consist of a superposition of emissions from many 
GRBs. In Section 5.1 we derive the radiation halo from a single powerful GRB "neutron bomb." 
Highly beamed GRBs will produce halos that are time averages of the emission spectra formed 
by a single GRB, weighted by the energy of events. Section 5.2 outlines energy deposition into 
the halo from neutron-decay protons. We also note that this process will make a very weak, long- 
lived afterglow of very energetic neutron-decay neutrinos. Section 5.3 describes multiwavelength 
prospects for detecting these halos. 



5.1. Radiation Halo from Neutron /3-Decay Electrons 

The neutron flux (neutrons cm~ 2 s _1 ) at location x in the case of a spherically symmetric 
explosion is simply 

®n(7n,t*;x) = exp(-x/c7 n t„) , (47) 

where i* is the time measured in the rest frame of GRB source and N(j n \ i*)^7n is the differential 
number of neutrons produced at time with Lorentz factors in the range between 7„ and 7 n + c?7n- 
The /3-decay electron and antineutrino each receive on average ~ 0.6 MeV from the decay (the 
neutron-proton mass difference is « 1.3 MeV). It is sufficiently accurate for the purposes here to 
let the proton and /3-decay electron each receive the same Lorentz factor as the neutron originally 
had. Thus the differential emissivity of either neutron-decay protons or electrons is simply 

(48) 

(see also Giovanoni and Kazanas (1990); Contopoulos and Kazanas (1995)). 

If we consider only the neutrons with j n » 10 4 which decay on timescales S> 10 7 s, then the 
GRB explosion and afterglow can be approximated as a 5-function in time. If we also approximate 
the neutron production spectrum as a single power law, then the neutron source spectrum can be 
represented by 

AT( 7n , U) = K nl - q 5{U - U)H[ ln ; 1, 7n , max ] , (49) 

where q is the spectral index of the neutron number spectrum, and 7 nimax ~ 10 10 is the maximum 
neutron Lorentz factor. Normalizing this spectrum to the total energy E n in neutrons, we have 

= ^ H^t) • (50) 

m n c 7n,max - 1 

Because q ~ 1, it is irrelevant whether the minimum value of 7„ is 1, as used here, or 10 4 . Without 
loss of generality, we let the explosion time = 0. Substituting equations (49) and (47) into 
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equation (48) gives 

K n 7~(9 +1 ) 

n(7,t*;x) = — - — — 5- i?[7;l,7n,max]exp(-a;/c7t n )(5(t* -x/c) . (51) 
cr„ 47rar 

We now consider the radiation signature of the neutron-decay electrons. Subsequent transport 
of the electrons can be neglected if an on-the-spot approximation is valid, which holds if the electron 
Larmor radius is much less than A n . This requires that the halo magnetic field B 3> m e c/(et n ) = 
6 x 10 -11 G. Faraday rotation measures of galaxy clusters and inferences from synchrotron radio 
halos indicate that cluster fields are ~ /J.G (see, e.g., the review by Eilek (1999)). It seems likely 
that galaxy halos would also be this strong. The solution to the electron continuity equation 



dt 



ID 

/°° P' An" 

H n{i, t>) , where t' = t- -^ 77 - . (53) 
I -at I 

Synchrotron radiation and Compton scattering of the cosmic microwave background radiation will 
dominate the energy losses of the electrons, although a more detailed treatment must treat Klein- 
Nishina effects which become important for electrons with 7 > (m e c 2 /fee ?{,{,) ~ 2 x 10 9 , where 
Tbb = 2.7 K. These loss rates can be written as —7 = i^o7 2 (see discussion following equation (23)). 
Substituting equation (51) into equation (53) and solving gives the result 

x, U) = ^ 2 exp(-x/ct„ 7 ') , (54) 

where 7' = [y^ 1 — z^o(t* — x/c)]' 1 and 

X <**<^ + ^ - 1 (7- 1 -7n,maJ- (55) 



Thomson and synchrotron losses can be treated on equal footing by rewriting equation (23) in 
terms of the photon emissivity 

n ph (e;x,t*) = ^-v e f ~ 1/2 e~ 3/2 n(< hz;x,t*) . (56) 



e 

The quantity e = en = B/B CT for synchrotron emission, and e = €{ for Thomson scattering, 
where is the characteristic dimensionless photon energy of the radiation field. The frequency 
= 4ca"Tiii/3, where the dimensionless field energy density m = ub for synchrotron emission and 
Ui = J °° cfcjejn so f t (ei) for Thomson scattering, where n so f t (ej) is the spectral density of the soft 
photon field. These equations are valid when <C 1. The restriction to the classical synchrotron 
regime always holds in this system, but the restriction to the Thomson regime may not apply, as 
already noted. 
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It is elementary to substitute equation (54) into equation (56) to obtain the photon emissivity 
ri p /j(e; x, t*) at location x and time t*. The spectrum observed at time t requires an integration 
over volume. Because the neutrons are flowing out at speeds very close to the speed of light, the 
expression t = t* + x(l — £t)/c accurately relates the explosion frame time and the observer time. 
Taking this relationship into account finally gives the synchrotron or Thomson spectrum observed 
at time t after a GRB explosion. It is 

vL v {t) (ergs s _1 ) = m e c 2 e 2 J dV h ph [e;x,U(t)] = ^^^^ \j\ 
r i rct/(2-n) 

x dfi dx [7(t)] 1_9 exp[-x/ct n j(t)] , (57) 

J-l ■/max{0, T 2f 7;j [t-^ 1 ( v /f-7^,Lax)]} 

where 

7W = {^-^-^(2-/x)]}- 1 . (58) 

The flux density S(v)(Jy) = 10 23 (vL l/ )/(4ird 2 L i') where v is the observing frequency. 

Fig. 8 shows calculations of the synchrotron and Thomson spectra emitted by neutron /3-decay 
electrons using equation (57). Here it assumed that 10 52 ergs in neutrons are emitted with a 
spectrum q = 1 up to a maximum Lorentz factor 7 nimax = 10 9 in Fig. 8a, and up to "f n ,ma,x = 10 11 
in Fig. 8b. In both calculations, we use a magnetic field B = 1 fiG and approximate the cosmic 
microwave background radiation as a <5-function soft photon source with dimensionless photon 
energy e = 4.6x 10~ 10 and energy density 4.1 x 10~ 13 ergs cm~ 3 . Note that although the synchrotron 
and Thomson spectra are plotted in the same graph, they are independently calculated. 

The peak luminosities reach ~ 10 36 ergs s _1 in Fig. 8a and ~ 10 34 ergs s _1 in Fig. 8b. 
The discrepancy with equation (46) implies that T ~ 0.1. This value is understood when one 
considers temporal smearing due to the finite energy-loss timescale, light travel-time effects and, 
most importantly, the contribution of late time (t 3> 7n,max^n) radiation. The bandwidth correction 
factor should also be considered. For a 1 //G field, the energy-loss timescale is ~ 7.7 x 10 20 7 _1 
s, which for 7 riimax = 10 9 is comparable to the 10 12 s neutron decay timescale. The temporal 
smearing due to light travel-time effects arising from emission produced on the far side of the 
explosion produces the high-energy features in the spectra observed at late times, particularly in 
Fig. 2b. In synchrotron and Thomson processes, the integrated luminosity decays oc t~ 1 at late 
times when most of the energy is injected in the form of high-energy electrons, as is the case here. 
Thus there is comparable energy radiated per decade of time at late times. A value of T ~ 0.1 in 
equation (46) is therefore reasonable. 

A notable feature in Fig. 8 is the appearance of sharp emission peaks at late times. These 
are the pileups that appear when electrons are injected with number indices harder than 2 and 
lose energy through synchrotron and Compton processes (Pinkau 1980). The synchrotron pileup 
features might be considerably broadened by magnetic-field gradients in the halos of galaxies. 
The situation regarding the pileup features in the Thomson peaks at ultra-high 7-ray energies is 
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more complicated and will require further study. Besides the Klein-Nishina effects that are not 
considered here and are crucially important in Fig. 2b, 7 rays with energies > 100 TeV ~ 10 28 Hz 
will materialize into e + -e~ pairs through 7-7 interactions with the cosmic microwave background 
radiation to form a pair halo surrounding the galaxy (Aharonian et al. 1994). Photons with energies 
> 10 TeV will not be observed due to pair-production attenuation on the diffuse infrared radiation 
field. The energy processed by this electromagnetic cascade will be transferred, in most cases, from 
the Thomson to the synchrotron components (consider, however, Kirk and Mastichiadis (1992)). 
Given detailed modeling and sensitive observations, the relative powers in the X-ray/soft 7-ray 
synchrotron component and the high-energy 7-component could, in principle, be used to infer the 
halo magnetic field. 

Effects of different magnetic field geometries and spatial variations of neutron injection due 
to beaming in GRBs should be considered in future work. Fig. 8 represents the simplest field 
geometry, but is representative of the integrated emission spectrum from radiation halos due to 
UHECR production in GRBs. 

5.2. Radiation Halos from Neutron-Decay Protons 

The neutron-decay protons carry three orders of magnitude more energy than the neutron- 
decay electrons, but this energy is also more difficult to extract. The proton Larmor radius is 
tl — IO710 B _1 (^G) kpc, so that much of the energy will be carried directly into intergalactic 
space when 7 nimax — 10 11 , even in the optimistic case of an extended (> 100 kpc) galaxy plasma 
halo with a mean magnetic field of ~ 1 ^uG. Such neutrons, together with the ultra-high energy 
protons and ions that diffusively escape from the blast wave, are of course postulated here to 
constitute the UHECRs. The Larmor timescale £l = nV c ~ 10 12 7ioB _1 (//G) s, so that a 10 19 eV 
proton might random walk for ~ 10 14 s before diffusively escaping from a 1 fiG, 100 kpc halo into 
intergalactic space. The timescale for energy loss through secondary production is (nhaio^ppc) 1 — 
10 15 /nhaio( cm_3 ) s ; where the mean halo particle density n^aio *C 1 cm -3 . Even though secondary 
production is very inefficient, it could however compete with the energy deposition by neutron 
/3-decay electrons if rihaio ^ 10~ 2 cm~ 3 , 7 riimax ^ 10 10 , and B > 1 ^G. Although the distinctive 
signature of secondary production is the 7r° — > 27 decay bump at ~ 70 MeV, it would be severely 
broadened due to the large Lorentz factors involved and would probably either not be detectable, 
or would form a low plateau to the diffuse galactic and intergalactic 7 radiation fields. 

Streaming instabilities excited by the outflowing neutron-decay protons could convert a large 
fraction of the available energy into long wavelength magnetohydrodynamic turbulence. Subse- 
quent cascades of the turbulence energy to shorter wavelengths could accelerate electrons through 
gyroresonant interactions with whistler and Alfven waves. Such processes have been invoked to 
explain the formation of diffuse radio halos in rich clusters (for a recent review, see articles in the 
collection edited by Bohringer et al. (1999)). It would be difficult, however, to distinguish neutron- 
decay halo nonthermal radio emission from electrons accelerated by cluster merger shocks (e.g., 
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Loeb and Waxman 2000). Searches for neutron-decay halos from field galaxies would therefore be 
more definitive than searching for such halos around galaxies within or near the peripheries of a 
galaxy cluster. 

5.3. Prospects for Detecting GRB Neutron-Decay Radiation Halos 

Three neutron-decay radiation halos are distinguished. In the Type-(3 halo, the power of the 
halo radiation field comes from /3-decay electrons. The previous section outlined in sufficient detail 
the principal radiative properties of a (3 halo. The most important uncertainty, besides the ever- 
present question of GRB source collimation, is the ratio of the magnetic-field energy density in 
the halo to that in the cosmic background radiation. The /3-decay electrons in a synchrotron (3 
halo place most of the radiated power in the synchrotron component. In contrast, microwave or 
ambient photons are Compton-scattered to ultra-high 7-ray energies to precipitate a pair shower 
in a Compton (3 halo. 

In the Type-p (for proton) halo, the power of the halo radiation field comes from /3-decay 
protons. A p-halo can be much brighter than a (3 halo, because it has a factor ~ 2000 more energy 
available, but the extraction and subsequent reradiation of this energy is far less easily quantified 
than for the (3 halo. Depending on the radiation transfer and environmental effects there are, as 
for the (3 halos, synchrotron p-halos and Compton p-halos. 

The third type of halo is the Typc-v (for neutrino) halo. The instantaneous neutrino energy 
spectra received at different times after the GRB can be obtained by following the approach of 
Section 5.1. The detection of a v halo is not technically feasible at present. 

5.3.1. Statistics of Neutron-Decay Halos 

Our starting point was the differential source rate density, equation (1). There we noted that 
in the power-law approximations for the Eq- and To-dependences of the differential rate density, 
one-half of the energy generated by the sources of GRBs comes from cosmic sources with apparent 
isotropic energy releases > 2.4 x 10 53 ergs. Fig. 3 shows that the neutron production efficiency 
increases monotonically with energy; therefore most of the neutron energy comes from GRBs with 
Eq > 3 x 10 53 ergs. These very energetic GRBs are, of course, much less frequent. Our study of 
GRB statistics (Bottcher and Dermer 2000) shows that in the universe on small (z < 0.1) scale, 
the rate density of GRB sources is ~ 3.6[E^2' 52 — (100)" ' 52 ] Gpc~ 3 yr _1 . Thus on average there 
are 0.43 GRB-type explosions per Gpc 3 per year with energy > 2 x 10 53 ergs. 

There are, speaking crudely, ~ Jl*/2 ^ ^K-^) — 0.52<I>* (L* galaxies)/Mpc 3 , so that the density 
of L* galaxies in the local universe is n^* ~ 2.3 x 10~ 3 Mpc~ 3 (compare Section 2.1). If all the 
mass of galaxies were wrapped up in L* galaxies, then each galaxy would see on average a GRB- 
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type explosion with energy > 2 x 10 53 ergs every ~ 5 Myrs, assuming that GRB explosions are 
uncollimated. Let t u represent the characteristic FWHM duration when the emission at frequency 
v from a neutron-decay halo reaches its peak luminosity Lq{u) = uL v . If t v <C 5 Myr, then the 
fraction of L* galaxies displaying emission at this level is ~ t u /5 Myrs. If, on the other hand, t v S> 5 
Myrs, then the galaxy will exhibit a superposition of the emissions from many GRB neutron-decay 
halos, with the total halo brightness reaching ~ Lo(v){t v /5 Myr). If the beaming fraction is ~ 1/500 
of the full sky (Frail et al. 2001), then GRBs will take place about once every ~ 10, 000 years in an 
L* galaxy, and will display neutron-decay halos at a level corresponding to the average GRB power 
multiplied by the neutron /3-decay production efficiency. Under these circumstances, the average 
bolometric power of an L* galaxy from (3 halos is at a level of 2.5 x 10 39 ergs s _1 x 4 x (2/1836) x 
1% xk c i = W 35 k c i ergs s" 1 , where the factor of 4 accounts for clean and dirty fireballs, and the 
factor k c i > 1 is a correction factor due to the enhancement of neutrino production in a clumpy 
medium. 



5.3.2. (3 Halos and p-Halos: Essential Features 

The essential features of a (3 halo produced by a single uncollimated GRB are given by the peak 
photon frequency v p k, the duration i^ur of peak luminosity L p k, and radial extent of the halo. For 
a synchrotron (3 halo, v p k ~ 3 x 10 20 .B(^G )7 2 (1 + z)~ l Hz, and idur ~ + z)j n ,maxt n = 3 x 10 5 7io 
yrs. Setting T ~ 0.1 in equation (46), the peak luminosity ~ 10 35 £ ; 547^ ) 1 ergs s -1 . The radial 
extent of the halo is ~ IOO710 kpc. 

A Compton /3-halo will be formed if the mean halo mag netic field (B) < 3(1 + zf pG. 
In this case, cosmic microwave background photons are Thomson scattered to energies ~ min(5 x 
10 3 , 5 x 10 4 7io)7io TeV. Many of these photons will materialize into electron-positron pairs through 
interactions with the cosmic diffuse background radiation field (Gould and Schreder 1967) to initiate 
an electromagnetic cascade that channels the radiant power into lower energy 7 rays and into a 
radially extended synchrotron component. The cascade ends when the photons penetrate the 
optical depth of the universe to 77 attenuation. This quantity is not well known, but Stecker and 
de Jager (1998) calculate that r 77 = 0.5 — 1 for ~ TeV photons from sources at z = 0.075 due to 
absorption by the diffuse intergalactic infrared radiation field. 

The p-halo will be brightest if the neutron-decay protons transfer and radiate their energy on 
a timescale shorter than the light-crossing time j n ,ma,xt n - Given the model-dependent uncertainty 
of the emergent photon spectrum from a p-halo, we approximate it with a vL v spectrum that has 
constant value L p between 10 6 Hz and 10 26 Hz. The uL u power radiated from the p-halo formed 
by a single strong GRB is, in this crude approximation for the spectrum, therefore at best 

/^(ergs s' 1 ) - ^± ~ 4 x 10 35 ^ , for 10 6 < u(Bz) < 10 26 , (59) 

ln(10 2U ) 710 710 

where we use T = 0.1 (compare eq.[46]). If the radiation is emitted in a narrow bandwidth, the p- 
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halo could be 2-3 orders of magnitude brighter. Thus there is emission in all observable wavebands 
at the level given by equation (59) during a period of ~ 10 13 7io s. The first-generation emission is 
distributed over a region of size 77 ~ IOO710 kpc, but the cascade radiation from the pair halo can 
occupy a much larger volume. 

5.3.3. Halo Detection: Observational Issues 

GRBs were first detected with soft 7-ray instruments (see Dermer et al. (1999)). To survey 
prospects for detecting neutron-decay halos, we begin at soft 7-ray energies and move to lower 
frequencies, returning at the end to the high-energy 7-ray domain. 

Soft "y-ray and X-ray Detection 

The sensitivity limit of a detector such as BATSE is ~ 0.2 x 100 keV cm' 2 s" 1 ~ 3 x 10~ 8 ergs 
cm~ 2 s _1 for a ~ 10-100 s observation, and that of OSSE is ~ 10 -11 ergs cm~ 2 s _1 for a two- 
week observation. Even with many orders of magnitudes improvement in sensitivity as provided by 
pointed instruments or position-sensitive technology, the detection of a neutron-decay halo is not 
easy with available X-ray detectors, much less 7-ray detectors. We estimate the limiting detection 
distance du m for a telescope with vF v sensitivity S = 10~ 15 S-is er g s over its nominal point- 
source observing time and bandpass. The peak luminosity of a neutron-decay (5 halo is given by 
equation (46), so that du m = ^ L/AnS = {E^/^iqS-i^) 1 / 2 Mpc. This criterion eliminates all 
7-ray instruments and all but the best X-ray detectors, such as Chandra 4 " with 5_i5 ~ 1. Within 
a few Mpc, the Milky Way and M31 are the closest L*-type galaxies. The rough odds are that 
a detectable halo could be observed from ~ 8710% of nearby L* galaxies if GRBs releases their 
energy isotropically, leaving only the two L* galaxy candidates if du m ~ 1 Mpc. A neutron-decay 
halo from M31 would cover a half-angle extent of O1/2 ~ 87° . Even for galaxies at ~ 10 Mpc, the 
challenge of background subtraction to reveal a cleaned X-ray image is severe, but would be assisted 
with model templates. It is worth recalling that beaming can increase the chance odds of sighting a 
galaxy that harbors a neutron-decay halo, but the halo itself would be at a proportionately smaller 
flux. 



4 Unfortunately, the fact that the neutron-decay halos are spread over a region greatly exceeding the extent of the 
galaxy makes them more difficult to detect, because they they will be harder to resolve from the diffuse background. 
The strength of instruments such as Chandra is that it focuses all photons from a point source onto one or a few 
pixels, so that the back-ground is greatly reduced (M. Bottcher, private communication, 2000). 
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Optical Detection 

In the spherical region that surrounds us to a depth of 100 Mpc, or within z = 0.022 for h = 0.65, 
there are, according to the earlier statistic, ~ 10 4 L* galaxies. At lOOdioo Mpc Mpc, the half-angular 
extent of a neutron halo is ~ 3.4(7io/dioo Mpc) arc minutes. At a sampling distance between ~ 10 
and 100 Mpc, there are therefore abundant candidates with galaxy disk sizes of ~ (2-20)7io arc 
seconds and a halo angular extent appropriate for an optical CCD. In the following, we sketch some 
basic considerations that enter optical halo detection. 

• The predicted /3-halo optical luminosity is « 10 35 ergs s _1 , but a neutron-decay halo could 
be as bright as 10 37 — 10 38 ergs s -1 if the parameters in the model are most optimistically 
tuned in favor of detecting a synchrotron p-halo. Compared to the typical L* galaxy optical 
luminosity of ~ 2 x 10 n LQ ~ 6 x 10 44 ergs s _1 , the halo luminosity provides a very weak 
flux. On the other hand, the emission is spread over a region that is far outside the optical 
radius of the galaxy. 

• The relative brightnesses of the central source and halo is ~ 6-9 orders of magnitude, or ~ 15- 
22 magnitudes. If the limiting magnitude is my = 25 for a good ground-based telescope, then 
a halo could only be seen for galaxies with my < 10. Noting that my ~ 5 for M31 implies 
that the limiting distance to detect a neutron-decay halo is ~ 10 Mpc for 2-3 meter class 
ground-based telescopes. In this case, the advantage of a halo that fills the CCD is lost, and 
the sensitivity of most large-aperture telescopes may not be good enough to detect the halo 
above background sources, the sky brightness and detector noise. 

• The limiting magnitude of the Hubble Space Telescope for point sources is my = 30. We 
could then potentially see neutron-decay halos to du m « 100 Mpc. Even at 100 Mpc, the 
halo subtends much of the CCD and the central bright source emission would have to be 
subtracted. For comparison, when subtracting central source flux from galaxy-disk flux in 
HST images, the contrast between the optical power of the AGN and that of the extended 
disk might have been ~ 10 2 -10 3 (this estimate is made by comparing optical luminosities of 
typical galaxies and QSOs, though the ratio could be even larger in studies where blazar light 
is subtracted from the host galaxy.) This still does not compare with the extreme contrast 
between the surface brightnesses of the optical disk of a galaxy and the surrounding diffuse 
halo. It seems that a blocking crystal for ground and space-based optical telescopes could be 
developed to eliminate the intense flux of the much brighter galaxy disks. The instruments on 
the Solar and Heliospheric Observatory probably achieve the greatest technical feat to detect 
faint objects in the field of a bright source (|m© — m sta rs| or \m & — m come t s | implies > 30 
orders of magnitude blockage of the Sun), but the detection of halos around distant galaxies 
will clearly pose different problems. 



Optical central-source luminosity is suppressed in certain classes of galaxies, most remark- 
ably, those that are likely to harbor active star formation. Here we are thinking of edge-on 
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starbursts (M82 or NGC 253-types) and dusty spirals, tidally-disturbed systems (e.g. Mrk 
421 and its satellite galaxies (Gorham et al. 2000)), and infrared luminous mergers such as 
Arp 220, Mrk 273, and other non-quasar members in Arp's atlas of peculiar galaxies. The 
search for neutron-decay halos also introduces a new avenue to examine the relative power of 
ULIGs (ultra-luminous IR galaxies) in stellar formation and black hole activity. The Infrared 
Space Observatory results on PHA/infrared line tracers of the starburst and AGN activity 
(Lutz et al. 1996) showed a separation of different galaxy types in a way that can be tested, 
because the strength of the neutron-decay halo is proportional to star-formation activity. The 
magnitude of either a (3 halo or a p-halo is, in this picture, directly proportional to the rate at 
which high-mass stars are formed, and is a basic assumption of the GRB statistics treatment 
of Bottcher and Dermer (2000). 

• AGNs and quasars introduce greater background subtraction problems, and pose the added 
difficulty of an interfering zodiacal light from high-latitude dust or gas that scatters the optical 
emission from the galaxy's AGN and stellar radiation fields. The existence of rather dense 
high latitude (~ 10-100 kpc) dust seems quite likely in an AGN environment due, for example, 
to tidal activity, disk winds, AGN radiation pressure on surrounding gas, and gravitational 
effects from distorted dark-matter halos and galaxy bars. Diffuse scattering plasma might 
also, unfortunately, be found in ULIGs for the same reasons. 

Technical considerations for detecting neutron-decay halos with optical telescopes will require 
an examination beyond the scope of this paper. A central insight is that even though point- 
source ("light-bucket") fluxes dim with source distance according to <fi oc d~ 2 in the Newtonian 
limit, the surface brightness of an optically thin source is constant (again, in the Newtonian limit). 
This effect has fundamental implications for observations against a source-confused and sky-limited 
background. 

Radio Detection 

The radio regime has the best vF v sensitivity, with S r ~ 10 9 Hz x 0.1 mJy ~ 10 -18 ergs cm~ 2 s _1 , 
combined with excellent angular resolution. Improved resolution (VLBI) must trade off with better 
limiting sensitivity ( VL A) , both of which additionally depend on observing frequency. For an opti- 
mistic radio halo power of ~ 10 35 ergs s _1 , the limiting sampling distance is only ~ 30(£ , 54/ / 7io) 1 / 2 
Mpc. Two effects determine the actual radio luminosity of a j3 halo. The first, as seen in Fig. 8, 
is that the radio luminosity is ~ 10 8 times dimmer than the peak nonthermal synchrotron power 
from a synchrotron (3 halo for our standard halo with a randomly oriented ~ 1 mean magnetic 
field. This reduction is partially offset by the fact that the synchrotron decay timescale from the 
radio-emitting electrons and positrons is larger by a factor of ~ 10 16 ' 5 /1.2 x 10 14 ~ 260 than the 
burst timescale (see Fig. 8) in the case of a 1 ^G halo. The net result is to reduce the sampling 
distance so that detecting the radio emission from a halo turns out again to be difficult. In the 
event of a very weak ((B) <C 0.1 fj,G) halo magnetic field, the reprocessing of the Compton power 
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into the synchrotron component could however improve radio detectability by moving u p k to lower 
frequencies (see Section 5.3.2). To take advantage of the good radio resolution (the size of the radio 
halo for sources at z ~ 1 (cz/Hq = 4600 Mpc) is on the order of a few arc-seconds) would require 
detection of (~ 10 30 ergs s 1 x 10 23 /10 57 cm 2 -10 9 Hz) < 10 7 //Jy fields spread over a surface area 
of this extent. This is not yet feasible. The tradeoff between angular extent and sensitivity will be 
helped if radio techniques can yield cleaned images that are ~ 0(°) in extent. Wherever the radio 
range proves to have the greatest capability (probably for galaxies at a few tens of Mpc), structure 
in the neutron-decay radio halo should be carefully sought. The advantage here is that to test the 
UHECR/GRB hypothesis, every L* galaxy should have a diffuse neutron-decay radio halo, whether 
or not GRB outflows are collimated. 

The very low-frequency (< 100 MHz) emission from neutron-decay halos persists around all 
L* galaxies, and forms part of the diffuse low-frequency radio background. Whether detection of 
such halos is technically feasible with new-generation radio arrays (e.g., the planned low- frequency 
array LOFAR) will require more study. 

High Energy j-ray Detection 

Returning now to the ultra-high energy gamma-ray regime, neither GLAST nor the ground-based 
air and water Cherenkov telescopes operating or in development (e.g., Whipple, Milagro, HESS, 
VERITAS) can be expected to detect a Compton (3 halo. Only under the most optimistic conditions 
of a highly luminous Compton p-halo at ~ 10 38 ergs s _1 is detection feasible. GLAST is ~ 50 times 
more sensitive than EGRET, which had a limiting sensitivity of « few x 10 -11 ergs cm~ 2 s _1 , as 
does Whipple. This gives a sampling distance of ~ 3OO-E54/710 kpc. With the seven- fold increase in 
limiting distance for GLAST, and with the improvement that will be achieved with the VERITAS 
array, there remains a chance of detecting highly luminous Compton p-halos from nearby galaxies. 

In summary, the search for direct synchrotron, and both direct and cascade 7 radiation from 
neutron-decay halos predicted by the external shock model are at limits that challenge current 
radio, optical, X-ray, and 7-ray detector technology. The predicted weakness of the /3-decay halos 
reflects the low neutron production efficiency in this model. The stronger neutrino and neutron 
production within the context of an internal shock scenario (Waxman and Bahcall 2000; Dai and Lu 
2001) would imply neutron decay halos ~ 100 times brighter. However, because of the uncertainty 
in the brightness of a p-halo, the detection of neutron-decay halos at the level of ~ 10 38 ergs s _1 
around L* galaxies would not discriminate between internal and external shock models. Radiation 
halos detected at the level of ~ 10 39 -10 40 ergs s _1 around L* galaxies would be inconsistent with an 
external shock scenario, but would require stong damping of the very energy neutron-decay proton 
energy in an internal shell model. The much stronger neutrino production in the internal shock 
model will also provide a clear discriminant between internal and external shock scenarios. 
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6. Cosmic Ray Production by GRBs 

As summarized in the Introduction, observations indicate that GRBs are associated with su- 
pernovae taking place in star-forming galaxies. Statistical analyses within the external shock model 
show that progenitor sources of GRBs inject a time- and space-averaged power > 2.5 x 10 39 ergs s _1 
into an L* galaxy. More general arguments also show that the power injected into the Milky Way 
by fireball transient sources with relativstic outflows is at the level of ~ 10 40:l=1 ergs s _1 (Dermer 
2000b). This is much smaller than the ~ 10 42 ergs s" 1 thought to be injected into the Galaxy by 
SNe of all types. It is also somewhat smaller than the galactic cosmic-ray luminosity of ~ 5 x 10 40 
ergs s _1 that is estimated to be required to power hadronic cosmic rays, depending on the assumed 
efficiency for accelerating hadronic cosmic rays. But the global cosmic-ray power estimate assumes 
that the locally observed cosmic-ray energy density is typical throughout the galaxy, and that tem- 
poral stochastic variations are not large. Both of these assumptions could be wrong (Hunter et 
al. 1997; Pohl and Esposito 1998). The derived FT power is sufficiently close to suggest that the 
progenitor sources of GRBs could power a significant fraction of the hadronic cosmic rays in the 
Galaxy. 

We additionally note several difficulties for the conventional view that cosmic rays are acceler- 
ated by supernovae in the galaxy (see Dermer (2000a) for further detail) : (i) Spectral signatures of 
the hadronic cosmic-ray component associated with ir° emission features, which carries ~ 30-100 
times as much energy as the leptonic cosmic-ray component, have not been detected unambiguously 
in the vicinity of SNRs (Esposito et al. 1996). (ii) The unidentified EGRET sources have not been 
firmly associated with SNRs (Romero et al. 1999), and several candidate SNRs are more likely 
to be associated with pulsars (e.g., Mirabal et al. 2001). (ii) TeV gamma rays are not detected 
from SNRs at the level expected from hadronic acceleration in SNR shocks (Buckley et al. 1998; 
Aharonian et al. 2001). (iii) The measured spectrum of the diffuse galactic 7-ray background is 
harder than predicted if the locally measured cosmic-ray proton spectrum is typical of other places 
in the Milky Way (Hunter et al. 1997). (iv) The origin of cosmic rays at and above the knee of the 
cosmic-ray spectum and the smooth transition at the knee are difficult to explain with a SN shock 
model (Lagage and Cesarsky 1983; Dermer 2001). 

We therefore suggest that cosmic rays originate from the subclass of SNe that are progenitors 
of GRBs. (Milgrom and Usov (1996) and Dar and Plaga (1999) have also suggested that cosmic 
rays might be accelerated by the sources of GRBs.) How frequent are these events compared to 
other types of SNe? We can evaluate this rate from the statistical study of Bdttcher and Dermer 

(2000) , or by modifying this study in view of the constant energy reservoir result of Frail et al. 

(2001) . 

The "supernova unit" 

^ r events . 1.3 x 10~ 5 . . 

SN « ' 101 °L e . B -io3 yr ' s -h^t- " (Gpc yr > <60> 

is defined in terms of the number of events of a given type per 10 10 Solar luminosities in the 
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blue band per century, recalling from Section 2.1 that J ga i,B — 7.6 x 10 16 /i7oL©,,B Gpc -3 , where 
hjo = 0.7h. For reference, note also that 1 SNu = 7 x 10~ 5 GEM, where the conversion factor to 
galactic events per Myr [GEM = # of events /(MW galaxy -10 6 yr); Wijers et al. (1998)] uses a 
Milky Way blue-band luminosity Lmw,b — 7r(15 kpc) 2 • 20L Q pc~ 2 = 1.4 x 10 10 L Qv b (Binney and 
Merrifield 1998; Scalo and Wheeler 2002). 

The local rate density of FTs calculated by Bottcher and Dermer (2000) when EV^ n = 10~ 4 
in equation (3) is h bd = 440 Gpc" 3 yr" 1 , implying a rate of ~ 80 GEM, or SNu(FT) 6 x 10" 3 
in supernova units. In galaxies of type Sbc-Sd, SNu(II) =■ 0.7(±0.35)/if and SNu(Ib/c) = 
0.14(±0.07)/if for Type II and Type Ib/c SNe, respectively (Cappellaro et al. 1997; Scalo and 
Wheeler 2002). Hence the ratio SNu(II)/SNu(FT) 120 ± 60, and SNu{Ib/c)/SNu{FT) 
24 ± 12. The FT rate density is, however, very sensitive to the number of low energy GRBs, and 
so could be « 30 times less frequent if E^ m = 0.1. 

If we accept the standard energy reservoir result of Frail et al. (2001), then the rate of the 
brightest GRBs will be > 500 times the observed rate. For the FT rate of tigrb(z = 0) = 3.6 
Gpc -3 yr -1 , which is the minimum rate necessary to fit the BATSE statistics within the external 
shock model, this implies a FT within the Milky Way ~ 4 x 10 _4 -10 -3 yr -1 , depending on the L* 
density. 5 Given the uncertainties associated with the rate and beaming estimations, we can expect 
that a GRB or FT will occur about once every 10 3 -10 4 yr in the Galaxy. Consequently ~ 1 in every 
10-100 SNe will display relativistic outflows including, in ~ 10-50% of the CtiSGS, Si GRB. GRBs 
are probably associated with the rarer Type Ib/c SNe that have lost their hydrogen envelopes. 
This prediction can be tested with a sample of many dozens of SNe Ib/c, where follow-on radio 
monitoring is used to identify relativistic outflows (Kulkarni et al. 1998; Weiler et al. 2000). A small 
fraction of SNe in our Galaxy should exhibit strong hadronic signatures associated with cosmic ray 
production. 

There are ~ 200-1000 black holes formed per Myr in the Milky Way if a black hole is formed 
by every GRB and FT. Over the 1.2 x 10 10 yr age of our galaxy, ~ 2 x 10 6 -10 7 black holes are 
thus formed. Gravitational deflection of the black holes off other stars and molecular clouds would 
increase the scale height of the older black holes to values exhibited by the older K and M stellar 
populations (e.g., Bahcall and Soneira 1980). A population of isolated black holes that accrete 
matter from the ISM will thus be formed by GRBs (Armitage and Natarajan 1999; Dermer 2000b). 
Some of the EGRET unidentified sources could originate from accreting, isolated black holes. 



Two measures of the density of galaxies like the Milky Way, which is assumed to be LMike, are used in this 
paper. The 2.3 x 10 6 Gpc" 3 fi gure refers to the number density of galaxies with luminosities between 0.5L* and 2L* 
(Section 5.3.1). The method of energy- weighting that was applied to the Schechter luminosity function in equation 
(5) implies an L* density of 4.6 x 10 6 Gpc -3 . 
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7. Summary and Conclusions 

We considered implications of the hypothesis that UHECRs are accelerated by the sources of 
GRBs in this paper. Here are the main points deduced from this study: 

1. The statistical study of GRBs in the external shock model (Bottcher and Dermer 2000) gives 
the rate density and emissivity of GRBs, and the number of GRBs with different Lorentz 
factor Tq and total apparent isotropic energy Eq. The event rate is very uncertain for weak 
GRBs, but most of the energy comes from the rare strong GRBs with Eq > 2 x 10 53 ergs. 
If GRBs are highly beamed, as suggested by Frail et al. (2001) and Panaitescu and Kumar 
(2001), then GRBs and fireball transients with Eq > 10 51 ergs could occur as frequently as 
once every 10 3 -10 4 yrs in the Galaxy. 

2. Formulae for the nonthermal synchrotron emission spectra produced by an external shock 
were used to derive the comoving nonthermal synchrotron photon spectra in the blast wave. 
This radiation provides target photons for the very energetic nonthermal particles. Neutrons, 
neutrinos, positrons, and pairs are formed as byproducts of 7-p and 7-ion interactions. 

3. Neutron and neutrino production spectra and light curves formed in photomeson interactions 
are readily derived in the external shock model. The neutrino flux from individual GRBs is 
far too weak to be detected by km 2 neutrino detectors, because most of the energy is carried 
by relatively few, very energetic neutrinos. GRBs might still contribute the major fraction 
of the diffuse neutrino background for neutrinos with energies > 10 16 eV. The energy carried 
away by neutrons from very energetic GRBs with E54 = Eq/W 54 ergs > 0.2 will exceed 1% 
of the total energy; reverse shock emission giving enhanced target photons could make the 
neutron-production efficiency even larger, though SSC processes might reduce it. 

4. Galaxies with GRB activity will be surrounded by neutron-decay halos formed by emissions 
from /^-electrons and neutron-decay protons. The halo size is ~ IOO710 kpc, where 710 is the 
Lorentz factor of the neutrons that carry most of the energy from the GRB. The shortest 
halo emission lifetime from a single GRB is ~ 3 x 10 5 yrs. The peak luminosity of a halo 
from a single smooth-profile GRB is ~ 10 35 i?54/7io ergs s _1 . Depending on the magnetic field 
strength in halos of galaxies, the neutron /3-decay electrons will produce a nonthermal syn- 
chrotron (3 halo with peak luminosities at optical/X-ray /soft 7-ray energies, and a Compton 
(3 halo at very high (GeV-TeV) 7-ray energies when the high-energy electrons Compton- 
scatter photons of the cosmic microwave radiation and induce a cascade. The p-halo formed 
by neutron-decay protons is more difficult to quantify and could be much brighter than the 
emission from a (3 halo. 

5. Because of sensitivity and imaging capabilities, prospects for detecting neutron /3-decay halo 
emission are best at optical and radio frequencies. For optimistic model parameters, it might 
also be technically feasible to detect these halos at X-ray and 7-ray energies. The subtraction 
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of the light from the bright central galaxy is a major obstacle to halo detection at optical 
frequencies. Approximately 8710% of L* galaxies should display a (3 halo from a single GRB 
near the peak of its luminosity output for unbeamed GRBs. If GRBs are highly beamed, 
then essentially all L* galaxies will be surrounded by such halos, but at a weaker average flux. 
The average bolometric neutron-decay (5 halo emission surrounding an L* galaxy IS ~ 

1Q 35 

ergs s^ 1 from smooth-profile fireball transients, but could be greater for GRBs occurring in 
inhomogeneous surroundings. 

6. The emissivity of the progenitor sources of GRBs in our Galaxy potentially provides a large 
fraction of the luminosity required to power the galactic cosmic rays. If cosmic rays are 
accelerated by the supernovae that collapse to form a GRB, about 1 out of ~ 20-100 supernova 
remnants in the Galaxy should have harbored a GRB, so that detection of a strong hadronic 
signature from a subset of supernovae should support this model for the origin of cosmic rays. 

In addition to the search for GRB neutron-decay radiation halos around galaxies, further 
progress on these problems will be achieved by searching for neutrino emission from GRBs. De- 
tection of high-energy neutrinos from smooth profile GRBs would probably rule out the external 
shock model for the prompt phase of GRBs. Another important study is to search for the sister 
classes of GRBs that do not trigger GRB detectors. Present GRB telescopes have strong triggering 
biases against dirty and clean fireballs transients predicted by the external shock model that can 
be remedied with appropriate slewing strategies and new detector designs (Dermer et al. 1999). 
Most important is to identify hadronic emission from SNRs through gamma-ray observations, and 
to determine if cosmic ray production is typical of all SNRs or, as suggested here, only a small 
subset of SNe that are associated with GRBs. 

I thank T. K. Gaisser, M. Bottcher, K. Weiler, S. Woosley, M. Leising, C. Kouveliotou, R. 
Berrington, K. Wood, T. Galama, and R. Schlickeiser for comments, and the referees for their 
reports. The work of CD is supported by the Office of Naval Research and the NASA Astrophysics 
Theory Program (DPR S-13756G). 

A. Synchrotron Radiation Limit on Maximum Proton Energy 

We derive an expression for the maximum proton energy that results from a competition 
between the acceleration energy-gain rate and the synchrotron energy-loss rate. For the acceleration 
rate, we assume that a particle cannot gain a large fraction of its energy at a rate more rapid than the 
gyrofrequency (Guilbert et al. 1983; Kirk 1991; de Jager et al. 1996; Rachen and Mesazaros 1998). 
This condition follows from general considerations for either first- or second-order acceleration. 
In either case, p' = fh(eB /2nm p c), where p' is the particle momentum in the comoving frame 
and /l < 1 is an unknown parameter. For first-order Fermi acceleration in the Bohm diffusion 
approximation, /l = 2tt/3+, where (3 + c is the upstream speed. For second-order Fermi acceleration, 
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/l = (87r/3)/3^(5i?/i?) 2 , where Pac is the Alfven speed of the scattering centers and (SB/B) 2 is 
the fractional energy density of resonant waves (Dermer 2001). 

The synchrotron radiation limit is obtained by balancing the acceleration rate with the syn- 
chrotron energy-loss rate given through equation (32). It easily follows that the maximum observ- 
able proton Lorentz factor is 

^ , r ( a) x io" , < Ar >::„ /4 g , . (A!) 



a T B ' [e B n(cm- 3 )] 1 /4 i + (4 T ) 3 / 16 [e^K 3 )]^ " 



The synchrotron radiation limit is implemented in the calculations with /l = 1. Thus, for example, 
the timescales of the highest energy protons shown in Figs, la and lb derive from this limit. Other 
effects that limit maximum particle energy through second-order Fermi acceleration in a GRB blast 
wave, such as available time and the requirement that the Larmor radius of the accelerated particle 
be less than the blast- wave width, are considered elsewhere (Dermer and Humi 2001). 



B. Analytic Estimate of Photomeson Production Efficiency 

The comoving differential photon energy density in a blast-wave geometry is 

e V(/)^(|)^,a„d e ^. (Bl) 



The comoving differential photon density is n'(e') = e'u'(e')/m e c 2 e' 2 , and the synchrotron emission 



is assumed to be isotropic in the comoving frame, so that n'(e',fi') = In'(e'). 



The time scale for significant energy loss by photohadronic reactions is given, following equation 
(29), by 

t' p ~\ ~ ! jH de' f dp'il - //)n'(e', m'Kt-t^") , (B2) 



we 



where a P1 ^{e") a S\p' - (1 - e^/i p e% a » 2 x 10~ 28 cm 2 , e A « 640, e" = j' p e'(l - //), and 
now use primes to refer particle Lorentz factors to their proper frame. The comoving time available 
to undergo hadronic reactions is t' = x/Tc. The quantity 

n _ t> „ a d 2 L e A t> [°° 3 

represents the efficiency to lose energy through photohadronic processes. 

In a fast cooling scenario, the photon energy at the peak of the vF v spectrum is 

2T 500 

€pk - (TT^) 657 ™" - TT~z fc ' e e(eBn ) 1/2 r| 00 . (B4) 

The threshold condition ^' p e' = e/\ requires acceleration of protons with observer-frame energies 
7 p = r 2 £A/[(l + z)e p k] when scattering photons with e « e p k- Because j p oc T -2 £ t 3 / 4 following 
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the prompt phase, demands on particle acceleration to scatter photons with e = e pk are more easily 
satisfied during the early episodes (t < t^) of a GRB. 



fe = fe pk {t)\—] a » ■ (B5) 



The vF v synchrotron flux can therefore be written as 

e 

We evaluate equation (B3) at if = t' d , and denote quantities at this time with hats. From equations 
(21) and (23) at e = e pk corresponding to 7 = 71 in a fast cooling scenario, 

& o ^ 3 k p e e E 6 j^r^n^ V e ^ergs^ 

^ fc - 4tt4 i 8vr CCJTj iV ^7™n - 4 47rd|t d " iU d 2 8 (l + z) cm 2 s ' 

where y = r eA/[(l + z)e pk 'j p ]. Hence f e = f tpk {e/e pk ) a " and a v = 1/2 for e < e pk , and a„ = 
(2 — p)/2 for e > e pk in the fast cooling limit (at smaller photon energies, though still above the 
synchrotron self-absorption frequency, a v = 4/3). From eq. (B3), 

f,~ — WW*" 1 (Bg) 

20 m e c 2 (l + zf{2 - a v )^x\e pk ' V ' 

Thus 

We (Dermer et al. 1999, a; Bottcher and Dermer 2000; Dermer 2000b) have previously shown 
that joint consideration of blast wave temporal and spectral characteristics of external shock emis- 
sion, levels of the diffuse background radiation, and triggering characteristics of 7-ray detectors 
cause burst telescopes to be biased in favor of the detection of GRBs with the prompt phase vF v 
peak frequency e pk = e^t, where e<± et is the photon energy of the telescope's largest effective area. 
For the BATSE telescope, €d e t — 0.1-1. Parameter sets A and B imply that e pk = 0.33/(1 + z) 
and 57 = 10~ 3 , and e pk = 0.42/(1 + z) and tj = 2x 10 -4 , respectively. The estimate for e pk com- 
pares favorably with the results of Figs, la and b. This estimates apply to protons with energy 
E p =■ rge A m p c 2 /[(l + z)e pk ] 9* 5 x 10 17 (r /300) 2 /[(l + z)(e pk /0.1) eV, which follows from the 
condition y = 1. 

From equation (B6) it is easy to show that in the fast cooling limit, fj oc El l2 Mr X)l \ and 

— 1/3 

E p at progressively higher proton energies. 

The spectral model used to fit the BATSE statistics is degenerate in the quantity tiqTq 
(Bottcher and Dermer 2000), and e pk is also degenerate in this quantity (eq. [B4]). But the pho- 
tomeson production efficiency 77 oc n^ 3 F^ 3 which, when the invariance uoFq is removed, indicates 
that 

fj oc ^/n^ . (B8) 

Based on the measured directional energy releases and the duration distribution of GRBs, To ~ 100- 
1000 and uniform CBM densities uq ~ 10~ 2 -10 2 cm~ 3 are implied. The lower range of the densities 
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are similar to values deduced from GRB afterglow fits (Panaitescu and Kumar 2001; Wijers and 
Galama 1999; Panaitescu and Kumar 2001a). Clouds or clumps in the external medium with 
densities no 3> 10 2 cm -3 are argued in this model to produce shorter duration spikes in GRB light 
curves, and so from equation (B8) could be more neutrino and neutron luminous. The conclusions 
of neutrino and neutron power calculated in this paper apply to smooth profile GRBs, and we 
predict no coincident neutrino fluxes from this type of GRB. 
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Table 1. Standard Parameter Sets 



Variable 


Quantity 


Parameter Set A a 


Parameter Set B a 


E 


Total Energy c (ergs) 


2 X 10 53 


2 X 10 53 


r 


Initial Lorcntz factor 


300 


300 


e e 


Electron Energy Transfer Parameter 


0.5 


0.1 


es 


Magnetic Field Parameter 


10- 4 


io- 1 


£max 


Maximum Particle Energy Parameter 


1 


1 


n 


Density of Surrounding Mcdium d (cm -3 ) 


10 2 


10 2 


V 


Nonthermal Particle Injection Index 


2.2 


2.2 


I 


Nonthermal Proton Energy Fraction 


0.5 


0.5 



a Paramctcr set giving good spectral fits to 7-ray luminous phase of GRBs 
b Parameter set giving good spectral fits to afterglow phase of GRBs 
c Apparcnt isotropic energy release 

d Surrounding medium is assumed to be uniform density 
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Fig. 1. — Ratio of comoving time to the energy- loss timescales for the processes p+j — > p+ir , n+n + 
(open circles) and proton synchrotron radiation (filled circles) as a function of the proton energy 
E measured by an observer. Curves are denoted by the base 10 logarithm of the observing time 
in seconds, (a) Parameters giving good fits to GRBs during the prompt phase with e e = 0.5 
and es = 10 -4 . Dotted lines show approximate expressions for the photomeson process given by 
equation (27). (b) Parameters giving good fits to GRBs during the afterglow phase with e e = 0.1 
and eg = 0.1. Other parameters are given in Table 1. 
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Fig. 2. — Instantaneous production spectra for neutrons (a) and neutrinos (b), labeled by the 
base 10 logarithm of the observing time in seconds. The prompt phase parameter set (A) is used 
in (a) and the afterglow parameter set (B) in (b). The thick solid and dotted curves give the 
time-integrated spectra for neutrons and neutrinos, respectively. 
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Fig. 3. — Production efficiencies for neutrons, denoted by "n," and neutrinos, denoted by V," as 
a function of apparent isotropic explosion energy £54 in units of 10 54 ergs. Prompt and afterglow 
parameter sets are denoted by labels (A) and (B), respectively. Neutron production efficiencies are 
also shown for the standard parameter sets, except now with Tq = 30 and Tq = 3000. Inset shows 
the dependence of the neutron production efficiency on the injection index p of the nonthermal 
protons for parameter sets (A) and (B), shown by the solid curve and dotted curve, respectively. 
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Fig. 4. — Diagram illustrating the temporal behavior of the rate at which neutrons and neutrinos 
with energies E = ma are produced by photomeson production in a GRB blast wave. The 
dimensionless time r = t/td, where td is the deceleration time. The lower triple-dotted dashed 
lines and the upper dot-dashed lines bound the production of particles due to the lack of available 
photons or the synchrotron radiation limit, respectively. The other regimes are a consequence of 
the spectral shape and evolution of the synchrotron soft photons, and the evolution of the energy 
in nonthermal protons. The interior tetragon defines the time during which the blast wave evolves 
in the strong cooling regime. 
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Fig. 5. — Calculations of neutron and neutrino production time profiles at 10 , 10 , 10 , and 
10 20 eV energies. Neutron production time profiles are given by the solid and dotted curves, and 
neutrino time profiles are given by the short-dashed and long-dashed curves for parameter sets (A) 
and (B), respectively. No neutrinos are produced at 10 20 eV for parameter set (B), as can be seen 
in Fig. 2b. 
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Fig. 6. — Calculations of the diffuse neutrino background flux produced by smooth-profile GRBs 
with parameters given by sets (A) (solid curve) and (B) (dotted curve), except that integrations 
are performed over Eq and Tq using equation (1) for the differential GRB emissivity. Inset shows 
the calculated GRB neutrino background for the smooth-profile GRBs. 
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Fig. 7. — Calculations of expected event rate of neutrinos from smooth-profile GRBs for a 1 km 2 
detector for parameter sets (A) and (B), shown by the solid and dotted curves, respectively. 




Fig. 8. — Calculations of spatially integrated synchrotron and Thomson spectra radiated by neutron 
/3-decay electrons emitted from a single GRB that radiates 10 52 ergs in neutrons. The neutron 
decay spectra is assumed to be represented by a power-law spectrum with dN/d^ n oc 7" 1 up to the 
maximum Lorentz factors of 7n,max = 10^ in panel (a) and 7n,max = 10 11 in panel (b). 



